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Abstract 

q ■ After a brief historical review of the emergence of QCD as the quantum field 

O ■ theory of strong interactions, the basic notions of colour and gauge invariance 
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are introduced leading to the QCD Lagrangian. The second lecture is devoted 
to perturbative QCD, from tree-level processes to higher-order corrections in 



^jjH , renormalized perturbation theory, including jet production in e + e annihila- 

tion, hadronic r decays and deep inelastic scattering. The final two lectures 



treat various aspects of QCD beyond perturbation theory. The main theme is 
effective field theories, from heavy quarks to the light quark sector where the 
spontaneously broken chiral symmetry of QCD plays a crucial role. 
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1 Introduction 

Why do we still study QCD after more than 30 years? 

• By decision of the Nobel Prize Committee in 2004 [1], QCD is the correct theory of the strong 
interactions. 

• The parameters of QCD, the coupling strength a s and the quark masses, need to be measured as 
precisely as possible. 

• Electroweak processes of hadrons necessarily involve the strong interactions. 

• In searches for new physics at present and future accelerators, the "QCD background" must be 
understood quantitatively. 

• Although QCD is under control for high-energy processes, many open questions remain in the 
nonperturbative domain (confinement, chiral symmetry breaking, hadronization, . . . ) 

• Last but not least, QCD is a fascinating part of modern physics. The lectures will therefore start 
with a brief historical review of the developments in particle physics in the sixties and early sev- 
enties of the last century. 

The following lectures were given to an audience of young experimental particle physicists. Although the 
lectures emphasize some of the theoretical aspects of QCD, the mathematical level was kept reasonably 
low. The first two lectures cover the basics of QCD, from the concepts of colour and gauge invariance to 
some applications of perturbative QCD. The last two lectures treat aspects of QCD beyond perturbation 
theory. The main theme is effective field theories, from heavy quarks to the light quark sector where the 
spontaneously broken chiral symmetry plays a crucial role. 

1.1 Historical background 

Particle physics in the early sixties of the last century was not in a very satisfactory state. Only for the 
electromagnetic interactions of leptons a full-fledged quantum field theory (QFT) was available. Quan- 
tum electrodynamics (QED) produced increasingly precise predictions that were confirmed experimen- 
tally. Nevertheless, the methodology of renormalization, an essential aspect of the perturbative treatment 
of QED, was not universally accepted. Even among the founding fathers of QFT, the dissatisfaction with 
"sweeping the infinities under the rug" was widespread. At the Solvay Conference of 1961, Feynman 
confessed [2] that he did not "subscribe to the philosophy of renormalization". 

What is the essence of this controversial procedure of renormalization that has turned out to be cru- 
cial for the shaping of QCD and of the Standard Model altogether? Specializing to QED for definiteness, 
three main steps are important. 

- Amplitudes A{j>i\ eo, mo; A) depend on the momenta pi of the particles involved, on the parame- 
ters eo , mo of the QED Lagrangian and on a cutoff A that cuts off the high-momentum modes of 
the theory. The cutoff is essential because A^pf, eo, mo; A) diverges for A — ► oo, rendering the 
result meaningless. 

- With the help of measurable quantities (cross sections, particle four-momenta) one defines physical 
parameters e(/j,) , m(/j,) that depend in general on an arbitrary renormalization scale \x. One then 
trades eo, mo for the physical e(//), m{fj) to a given order in perturbation theory. 

- The limit lim A(pf, eo(e, m, A), mo(e, m, A); A) = A(pf, e(u), m(//)) is now finite and unam- 

A^oo 

biguous for the chosen definitions of e(//), m(/x). 

Based on this procedure, the agreement between theory and experiment was steadily improving. 

For the weak interactions, the Fermi theory (in the V — A version) was quite successful for weak 
decays, but 



• higher-order corrections were not calculable; 

• for scattering processes the theory became inconsistent for energies E > 300 GeV (unitarity 
problem). 

The rescue came at the end of the sixties in the form of the electroweak gauge theory of the Standard 
Model. 

Of all the fundamental interactions, the strong interactions were in the most deplorable state. 
Although the rapidly increasing number of hadrons could be classified successfully by the quark model of 
Gell-Mann and Zweig [3], the dynamics behind the quark model was a complete mystery. A perturbative 
treatment was clearly hopeless and the conviction gained ground that QFT might not be adequate for the 
strong interactions. 

This conviction was spelled out explicitly by the proponents of the bootstrap philosophy (Chew et 
al.). Under the banner of nuclear democracy, all hadrons were declared to be equal. Instead of looking for 
more fundamental constituents of hadrons, the S-matrix for strong processes was investigated directly 
without invoking any quantum field theory. Although the expectations were high, nuclear democracy 
shared the fate of the student movement of the late sixties: the promises could not be fulfilled. 

A less radical approach assumed that QFT could still be useful as a kind of toy model. The 
main proponent of this approach was Gell-Mann who suggested to abstract algebraic relations from a 
Lagrangian field theory model but then throw away the model ("French cuisine program" [4]). The 
usefulness of this approach had been demonstrated by Gell-Mann himself: current algebra and the quark 
model were impressive examples. Until the early seventies, Gell-Mann took his program seriously in 
declaring the quarks to be purely mathematical entities without any physical reality, a view shared by 
many particle physicists of the time. 

The decisive clue came from experiment. Started by the MIT-SLAC collaboration at the end of the 
sixties, deep inelastic scattering of leptons on nucleons and nuclei produced unexpected results. Whereas 
at low energies the cross sections were characterized by baryon resonance production, the behaviour 
at large energies and momentum transfer was surprisingly simple: the nucleons seemed to consist of 
noninteracting partons (Feynman). Obvious candidates for the partons were the quarks but this idea led 
to a seeming paradox. How could the quarks be quasi-free at high energies and yet be permanently bound 
in hadrons, a low-energy manifestation? 

That the strength of an interaction could be energy dependent was not really new to theorists. 
In QED, the vacuum acts like a polarisable medium leading to the phenomenon of charge screening. 
However, contrary to what the deep inelastic experiments seemed to suggest for the strong interactions, 
the effective charge in QED increases with energy: QED is ultraviolet unstable. 

To understand the phenomenon of an energy dependent interaction, we consider the dimensionless 

ratio of cross sections 

a(e + + e - — > hadrons) 

e+e~ / _j_ J Z Z J ZZ ■ v-U 

Beyond the leading-order value Rq (cf. Sec. 2.1), one finds to lowest order in the strong coupling 
constants g s : 

R e+ e- = R (l + |L) • ( 2 ) 

The general form to any order in g s (neglecting quark masses) is 

R e+e - = R e+e -(E,n,g s (n)) (3) 

where E is the center-of-mass energy and \x is the renormalization scale. Since R e + e - is a measurable 
quantity, it must be independent of the arbitrary scale fx: 

d ( d d \ 

fi—R e+e -(E, f i,g s ( y n)) = — ► U— + (3(g s )—\R e+e _ =0, (4) 



with the beta function 

a, x dg a (ji) 

P{9s) = n—j • (5) 

Dimensional analysis tells us that the dimensionless ratio R e + e - must be of the form 

E 
B e + e -{E,/j,,g s (ij,)) = /(— ,5s(m)) • (6) 

The seemingly uninteresting dependence on /x can therefore be traded for the dependence on energy or 
on the dimensionless ratio z = E//J, : 

+ /%,)*- )f(z,9s(fj)) = 0. (7) 

siogz ^ s / 

The general solution of this renormalization group equation is 

f{z,g s {n)) = f(g!(z,g 3 )) , (8) 

i.e., a function of a single variable, the energy dependent (running) coupling constant ~gl{z, g s ) satisfying 

dgl 



dlogz 



am (9) 



with the boundary condition g s (l,g s ) = g s . For any gauge coupling, the leading one-loop result for the 
(3 function is 

A. 3 

(4tt) 2 
implying 



ft*) = -TTZW^ 3 (10) 



gf{E/n,g s (n)) = R g ' (/i) • (ID 



^0 r 2,,,M ,.,2 i.;2 

(4vr) 2 - 



1 + TT^STs (M) log £7m" 



Expanding the denominator, we observe that the renormalization group equation has allowed us to sum 
the leading logs (gg(fi) logE 2 / ' jx 2 ) of all orders in perturbation theory. Even more importantly, the 
energy dependence of the running coupling constant is determined by the sign of j3q in Eq. (10): 



fa < 0: lim g(E) = infrared stable (QED) 

E — *0 

j3 > 0: lim g(E) = ultraviolet stable (QCD) 

E^QO 



R e+e - - Bo ( 1 + 4^ + 0{g A s {E)) ) (12) 



For the cross section ratio R e + e - we get finally 

J ~"47T 2 ' WV ' 7 '- 

in terms of g 2 s (E) =gf{E/n,g s {n)). 

The crucial question in the early seventies was therefore whether QFT was compatible with ultra- 
violet stability (asymptotic freedom)? The majority view was expressed in a paper by Zee [5]: "... we 
conjecture that there are no asymptotically free quantum field theories in four dimensions." While Cole- 
man and Gross set out to prove that conjecture their graduate students Politzer and Wilczek (together 
with Gross) tried to close a loophole: the P function for nonabelian gauge theories (Yang-Mills theories) 
was still unpublished and probably unknown to everybody except t'Hooft. In the spring of 1973, the 
Nobel prize winning work of Politzer and Gross and Wilczek [6] demonstrated that Yang-Mills theories 
are indeed asymptotically free. 



The crucial difference between QED and QCD is that photons are electrically neutral whereas the 
gluons as carriers of the strong interactions are coloured. Further physical insight can be obtained by 
taking up an analogy with the electrodynamics of continuous media [7]. Because of Lorentz invariance, 
the vacuum of a relativistic QFT is characterized by 

en = 1 (13) 

for the product of permittivity e and permeability \x. In QED charge screening implies e > 1 so that the 
vacuum of QED acts like a diamagnet (// < 1). In QCD the colour charge screening of quarks (e > 1) is 
overcompensated by gluons (spin 1) acting as permanent colour dipoles (/j, > 1). Because 

/3 QCD = i (11JV C - 2N F ) , (14) 

the QCD vacuum is a (colour) paramagnet for Np < llN c /2 < 17 quark flavours (for N c = 3). 
Because of the general relation (13) this can also be interpreted as anti-screening (e < 1). 

The existence of three colours was already widely accepted at that time. Gell-Mann and collabora- 
tors had been investigating a model of coloured quarks interacting via a singlet gluon (not asymptotically 
free). In a contribution of Fritzsch and Gell-Mann in the Proceedings of the High Energy Conference in 
Chicago in 1972 [8] one finds the probably first reference to nonabelian gluons: "Now the interesting 
question has been raised lately whether we should regard the gluons as well as the quarks as being non- 
singlets with respect to colour (J. Wess, private communication to B. Zumino)." Although Gell-Mann 
is generally credited for the name QCD, the first published occurrence of QCD is much less known 
(cf., e.g., Refs. [1]). My own investigation of the early literature has produced a footnote in a paper of 
Fritzsch, Gell-Mann and Minkowski in 1975 [9] suggesting "A good name for this theory is quantum 
chromody namics ." 

1.2 Colour 

Already before the arrival of QCD, there were a number of indications for the colour degree of freedom. 

• Triality problem 

In the original quark model often called the naive quark model, the three quarks u, d, s give 
rise to mesonic bound states of the form qq. All of the nine expected bound states had already 
been observed suggesting an attractive force between all quarks and antiquarks. The baryons fit 
nicely into qqq bound states. If the strong force is purely attractive why do antiquarks not bind 
to baryons? The resulting objects of the form qqqq would have fractional charge and have never 
been observed. Introducing three colours for each quark and antiquark allows for 9 x 9 = 81 
combinations of qq only nine of which had been found. The remaining 72 combinations are not 
bound states invalidating the previous argument. 

• Spin-statistics problem 
Consider the state 

|A ++ (S z = 3/2))~|uT u| «T>- (15) 

Since the spin-flavour content is completely symmetric, Fermi statistics for quarks seems to require 
an antisymmetric spatial wave function. On the other hand, for every reasonable potential the 
ground state is symmetric in the space variables. Colour solves this problem because the state 
(15) is totally antisymmetric in the colour indices respecting the generalized Pauli principle with a 
spatially symmetric wave function. 

• Renormalizability of the Standard Model 

With the usual charges of quarks and leptons, the Standard Model is a consistent gauge invariant 
QFT only if there are three species of quarks in order to cancel the so-called gauge anomalies. 



7T° — ► 2 7 decay 

The by far dominant contribution to the decay amplitude is due to the chiral anomaly (exact for 
massless quarks). The observed rate can only be understood if there are again three species of 
quarks. 

Momentum balance in deep inelastic scattering 

The momentum sum rule indicates that only about 50 % of the nucleon momentum is carried by 
valence quarks. The remainder is mainly carried by gluons. 

Quark counting 

After the advent of QCD, many more direct confirmations of the colour degree of freedom were 
obtained. One of the first confirmations was provided by the total cross section a(e + + e~ — ► 
hadrons) already discussed in the previous subsection. 

Hadronic r decays 

As we shall discuss in the next lecture, hadronic r decays not only give clear evidence for N c = 3 
but they also provide an excellent opportunity for extracting a s = g 2 /^- 



T(r -^v T + hadrons) 

- i itr = ^r, n (to) 

e , V- , d e r(r- -► v T e~V^ 




N c {\V ud \ 2 + \V us \ 2 )(l + 0(a s )) . 



w u 



Fig. 1: Feynman diagram for r — ► v T + X. 



1.3 Gauge invariance 

Gauge invariance is a main ingredient not only of QCD but of the Standard Model as a whole. We start 
with the Lagrangian for a single free Dirac fermion: 

Cq = ip(x) i$ ip(x) — mip(x) ip(x) , (ji := 7 M a M . (17) 

This Lagrangian and the resulting Dirac equation are invariant under a phase transformation (global 

u(i)) 

ip(x) — ► ip'(x) = e- l Q £ tp{x) (18) 

with Qe an arbitrary real constant. One may now pose the question whether the phase in the trans- 
formation law (18) must really be the same here and "behind the moon", as is the case in (18) with a 
space-time independent phase Qe. Instead of experimenting behind the moon, we replace the constant e 
with an arbitrary real function e(x) and see what happens. As is easily checked, the mass term in (17) 
remains invariant but not the kinetic term because 

d^{x) — ► e-'QsO*) {dp - iQd^s{x)) Y>(x) . (19) 

The conclusion is that the phase must indeed be the same here and behind the moon for the theory to be 
invariant under transformations of the form (18). 

However, there is a well-known procedure for enforcing local invariance, i.e. invariance for a 
completely arbitrary space-time dependent phase e(x). We enlarge the theory by introducing a spin-1 
vector field A^ that has precisely the right transformation property to cancel the obnoxious piece in (19) 
with d^e(x). The idea is to replace the ordinary derivative d^ by a covariant derivative D^: 

D^{x) = (d^ + iQA^x))^), (20) 



with 

A„(x) — ^(aO = A^x) + d^e(x) . (21) 

It is easy to check that D^ip transforms covariantly, 

D^(x) — {p,^)'{x) = e-'Q^D^x) , (22) 

so that the enlarged Lagrangian 

£ = lp(x) (ilp-m) ip(x) =C - QA^{x)$(x) yty(x) (23) 

is invariant under local 17(1) transformations (gauge invariance). 

The most important feature of this exercise is that the requirement of gauge invariance has gener- 
ated an interaction between the fermion field ip and the gauge field A^. Introducing a kinetic term for 
the gauge field to promote it to a propagating quantum field, the full Lagrangian 

£=^(ilp-m)ip- \f^ v F^ (24) 

is still gauge invariant because the field strength tensor F^ u = d^A v — d u A^ is automatically gauge 
invariant. Setting Q = — e for the electron field, we have "deduced" QED from the free electron theory 
and the requirement of gauge invariance. For completeness, we take note that a mass term of the form 
MlA^A^ is forbidden by gauge invariance implying massless photons. 

1.4 SU(3) C and the QCD Lagrangian 

As we have seen, quarks come in three colours. Suppressing all space-time dependence, the free quark 
Lagrangian for a single flavour has the form 

3 

A) = ^ Qi {i0 - m q ) qi . (25) 

Assuming the three quarks with different colours to have the same mass m q (different flavours have 
different masses, of course), we can ask for the global invariances of Cq. All transformations that leave 
Cq invariant are of the form 

qi^q'i = U ijqj , utf = U J <U=t (26) 

with arbitrary unitary matrices Uij. Splitting off a common phase transformation qi — ► e~ %£ qi treated 
previously and generating electromagnetic interactions that we know to be colour blind, we are left with 
the special unitary group SU(3) comprising all three-dimensional unitary matrices with unit determinant. 

In contrast to the £7(1) case treated before, we now have eight independent transformations in 
SU(3) (an 8-parameter Lie group). Therefore, continuing in the same spirit as before, it is not just a 
question of demanding gauge invariance but we also have to find out which part of SU(3) should be 
gauged. With hindsight, the following two criteria lead to a unique solution. 

i. The three colours are not like three arbitrary electric charges but are instead intimately connected 
through gauge transformations. This requires the quarks to be in an irreducible three-dimensional 
representation leaving only two possibilities: either all of SU(S) or one of the SU(2) subgroups 
must be gauged. 

ii. Quarks and antiquarks transform differently under gauge transformations (3 7^ 3*). This closes 
the case and implies that all SU(3) transformations must be gauged. 



The above requirements guarantee that both qq and qqq contain colour singlets = hadrons, 

3* 03 = 1© 8, 3©3©3 = lffi8ffi8©10, (27) 

but neither qq nor qqqq or other exotic combinations. 

Every three-dimensional unitary matrix with det U = 1 can be written as 

8 A 
U(e a )=exp{-iJ2e a Y} ( 28 ) 

a=l 

with eight parameters s a and with eight traceless hermitian Gell-Mann matrices A a . Their commutation 
relations define the Lie algebra of SU(3): 

[X a , X b ] = 2i f abc Xc , (29) 

with real, totally antisymmetric structure constants f abc . 

The gauge principle demands invariance of the theory for arbitrary space-time dependent functions 
e a (x). Instead of a single gauge field A^, we now need eight vector fields G^(x) entering the covariant 
derivative 

(0"g)i= (d»6 ij +ig 8 f^G£^jq j =:{(&> + ig 8 G*)q} i . (30) 

The real coupling constant g s measures the strength of the quark gluon interaction just as the charge 
Q is a measure of the electromagnetic interaction. Using the convenient matrix notation (summation 
convention implied) 

,A a 



G% ■= GZ^f- , (31) 



the covariant derivative now transforms as 



G M — G'„ = U(e)G^U\e) + - (0„l7(e)) U\e) (32) 

9s 

in order for (D^q)i to transform like the qi themselves (co variance requirement). 

Because of the nonabelian character of SU(3), the transformation laws are more complicated than 
in the electromagnetic case. The differences can already be seen in the infinitesimal transformations of 
the gluon fields G% following from (32): 

G^ — G% = G% + -d^s a + f abc e b G» + 0(e 2 ) . (33) 

9s 

In order to have propagating gluon fields, we need an analogue of the electromagnetic field strength 
tensor F^ v . The simplest approach is to calculate the commutator of two covariant derivatives: 

[D^ D v ] = [d^ + i gsGf,, d u + i g s G u ] =: i g s G^ u . (34) 

The nonabelian field strength tensor G^ v = G£ u — has the explicit form 

G> xv = g^c u -d u G'" + ig s [G >1 ,G u } (35) 

G? = d^ a -d^-g s f abc G^ c 



and it transforms covariantly under 577(3) gauge transformations: 

GV — G 1 ^ = U{e)G^U\e) . (36) 

The gauge invariant colour trace 

tr(GVGH = ^G^ (37) 

has the right structure for a gluon kinetic term leading immediately to the SU(3) C invariant QCD La- 
grangian for / = 1, . . . , Np quark flavours: 

1 N F 

£qcd = --XxiG^GT) + ^q f {0 ~ m/l c ) q f . (38) 

/=i 

As in the £7(1) case, gauge invariance requires massless gluons. Writing out the Lagrangian (38) in 
detail, one finds three types of vertices instead of a single one for QED: 

qqG GGG GGGG 




(TCCCCK TtnrcG-Q 



CV. /to 

c\ (.6 



Fig. 2: Basic vertices of QCD. 

In addition to the quark masses, QCD has a single parameter describing the strength of the strong inter- 
actions, the strong coupling constant g s (a s = gl/Air). 



Experimental group theory 



Can experimentalists determine more than a single coupling strength a s in a strong process? All the 
information is contained in the vertices: in addition to g s , the vertices also contain the two matrices 

(*o )ij = o (^ a )*i ' (*a )bc = -ifabc , (39) 

defining the fundamental and adjoint representations of (the Lie algebra of) 577(3): 

[ta,t b ] = i fabcU ■ (40) 

Let us pretend for a moment that we don't know that there are three colours and eight gluons. For a 
general (compact Lie) group of symmetry transformations, the vertices are again determined by quark 
and gluon representation matrices t a ,t a . The combinations that actually appear in measurable quantities 
are the following traces and sums: 

tv(t R t R )=T R S ab , Y.^M t a)ik = C R 5 ik (R = F,A), (41) 



with 



Tr\ Dynkin index for the representation R; 

Cr: (quadratic) Casimir for R. 



For a (^-dimensional representation, one derives from the definitions (41) the general relation 

d R C R = n G T R (42) 

where uq is the number of independent parameters of G. Restricting the discussion to SU(n), the two 
cases of interest are 

R = A (adjoint representation): 

dA = nG — ► Ca = Ta = n for SU(n); 

R = F (fundamental representation of SU(n)): 

- 2 - 1 r„ - 1/2 — ► C F 

4 



d,F = n, uq = n — 1, Tp 
and for the special case of 577(3): 



n 



1 



C F 



Ca = Ta 



2n 
N r -- 



The independent quantities that can be measured are Cf and Ca- A combined jet analysis in e + e 
annihilation at LEP found [10] 

C> = 1.30±0.01(stat)±0.09(sys), C A = 2.89 ± 0.03(stat) ± 0.21(sys) (43) 

in manifest agreement with 51/(3). 

Feynman diagrams are constructed with the vertices and propagators of the QCD Lagrangian (38). 
The problem here is the same as in QED: due to the gauge invariance of (38), the gluon propagator does 
not exist. At least for perturbation theory, the inescapable consequence is that gauge invariance must be 
broken in the Lagrangian! Or, in a more euphemistic manner of speaking, the gauge must be fixed. In the 
simplest and widely used version (co variant gauge with real parameter £) the Lagrangian (38) is replaced 

by 

e 



-QCD 



£qcd - 7T {dfi,G%) + £ g host • 



The gluon propagator now exists (£ = 1: Feynman gauge): 



AM 



lab- 



<r + (r 1 - 1 



k^k 1 



c=ij5 -tg- 

= o ab 



fXU 



(44) 



(45) 



k 2 + ie \ " k 2 J k 2 + ie 

The additional ghost Lagrangian £ g host repairs the damage done by gauge fixing: although Green func- 
tions are now gauge dependent, observable S-matrix elements are still gauge invariant and therefore 
independent of £ (Feynman, Faddeev, Popov, BRST, . . . ). 



2 Perturbative QCD 
2.1 QCD at tree level 

The calculation of tree amplitudes in QCD is straightforward but 

• to compare theory with experiment, we must have hadrons rather than quarks and gluons in the 
initial and final states; 

• amplitudes and cross sections are in general infrared divergent for massless gluons. 

The general recipe is to consider infrared safe quantities, the more inclusive the better. A good example 



is once more e + e 



hadrons. 




Fig. 3: e+e 



+ P - 



hadrons. 



hadrons 



in 



The sum over all hadronic final states can be expressed in terms of the imaginary (absorptive) part of the 
two-point function of electromagnetic currents (photonic case), the hadronic vacuum polarization: 



UZ(q) = i / d^xe^ • x (0\ TJ^ m (x)^ m (0)|0) = {-g^q 2 + gV) IW? 2 ) . 



(46) 




Fig. 4: Hadronic vacuum polarization. 



The sum over all intermediate states can be performed either with quarks and gluons or with hadrons. 
Since there are no massless hadrons, the hadronic vacuum polarization is infrared safe. 

To lowest order in QCD, the amplitude for 



e e 



7 *(Z* 



qq 



(47) 



is in fact independent of the strong coupling constant g s . Except for the charges, masses and multiplicities 
of quarks, it is the same calculation as for e + e~ — ► /J, + fJ,~ in QED. Therefore, for quarks with given 
flavour / and colour i the amplitude is (neglecting m e , m M , m q ) 



A{e + e~ -► q)q)) = ^A{e + e 



H 1 H 



(48) 



Quarks and antiquarks with different colour and flavour are in principle distinguishable so that the total 
hadronic cross section, normalized to a(e + e~ — ► /i + ^~), is 



Rp 



a(e + e~ 



hadrons) 



i,f f 



a(e + e — > /i + /x ) 
As shown in Fig. 5, this is a good approximation to the experimental data between quark thresholds. 



(49) 




Fig. 5: Experimental data for R e + e - taken from Ref. [11]. 



1 I 



A similar result is obtained for the hadronic width of the Z: 

R z = V{Z -► hadrons)/r(Z -► e+e") = 7V C (1 + few) ^( w / + a ))K v l + °: 

/ 

where vp, ap are the (axial-) vector couplings for Z — ► FF. 



2), 



(50) 



e ' e 



jets 



At high energies, the two-jet structure from e + e~ — ► gg dominates, being the only process at O(ag). At 
0(a s ) and omitting Z exchange, we have in addition gluon bremsstrahlung off quarks giving rise to a 
three-jet structure: 

e + {qi)e~{q 2 ) -► q{pi)q(p2)G(p 3 ) . (51) 





Fig. 6: Leading-order diagrams for three-jet production. 



The calculation is again identical to QED bremsstrahlung except for a factor (sum over all final states in 
the rate) 



^ tr(iftf) =T F '^2 &aa = T F nsu(3) = d F C F = 3C F = 4 . 



(52) 



With the kinematics specified by 

s = (qi + q2) 2 , (pi+Pj) 2 = (qi + 92 -Pkf =■ s(i -x k ) 



(53) 



^l + ^2 + a^3 



CMS : Xj = 2Ej 



the double differential cross section (for massless quarks) is found to be 

d 2 a 2a s cr x\ + x\ ... 47ra 2 



cixi dx2 



37T (1 — xi)(l — x 2 ) 



with 



00 






(54) 



The problem with this cross section is that it diverges for x\ — ► 1 (i = 1, 2). This infrared divergence is 
due the singular behaviour of the quark propagator and it happens even for massive quarks: 



(P2+P3) 



m n 



2^2 -P3 = s(l -xi) . 



(55) 



m g > 



m„ 







xi — > 1 only possible for P3 — > (soft gluon singularity); 
x\ = 1 also possible for P3 1 1 P2 (collinear singularity) . 



To understand the origin of infrared divergences, we first take the viewpoint of an experimentalist mea- 
suring three-jet events where the jets stand for the quarks and the gluon in the final state. 

• Depending on the detector resolution, a quark and a soft gluon cannot be distinguished from a 
single quark. In that case, the event will be counted as a two-jet event. 

• Two collinear massless particles can never be resolved: they always stay together. 
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From the viewpoint of a theorist, we recall that perturbation theory is built on the assumption that par- 
ticles do not interact when they are sufficiently far apart. This assumption is not really satisfied for 
massless quanta like photons or gluons that give rise to long-range forces. In other words, an electron (a 
quark) can never be separated from its cloud of soft photons (gluons). 

The practitioner's solution of the infrared problem is well understood: 

• One must define criteria to distinguish between (in the present case) two- and three-jet events (jet 
algorithms). 

• Virtual gluon (loop) corrections for the process e + e~ — ► qq must be included. 



2.2 Higher-order corrections and renormalization 

The loop corrections of 0(a s ) for e + e~ — ► qq are calculated from the Feynman diagrams below. 





e+ q 

Fig. 7: One-loop diagrams for e + e~ — ► qq. 

The resulting amplitudes are both infrared and ultraviolet divergent. 

In contrast to infrared divergences, ultraviolet divergences are due to the high-momentum compo- 
nents of the particles in loops. Common sense tells us that those components cannot influence physics at 
low energies. If this were the case we would have to give up all hopes of being able to make predictions 
at presently accessible energies. 

The recipe to handle ultraviolet divergences is also well understood. One first has to choose a 
method to cut off the high-momentum components. There are infinitely many ways to do that so the 
question is legitimate whether the final amplitudes will depend on that procedure rendering the result 
completely arbitrary. The answer is that the cutoff procedure (regularization) must always be accom- 
panied by renormalization. Before choosing a suitable regularization procedure let us therefore try to 
understand the idea of renormalization, using the most naive regularization method. 

To simplify matters as much as possible, we consider the elastic scattering of two particles in 
massless scalar <fi 4 theory (£j nt ~ A</> 4 ): 

H - H , (56) 

with scattering amplitude A(s, t) in terms of the usual Mandelstam variables. We now define what we 
mean by the physical (renormalized) coupling constant. The definition should be applicable at every 
order of perturbation theory and it should coincide with the constant A in the Lagrangian at tree level. A 
possible definition in scalar 4> 4 theory is 

Ar(/i) := A(s = ~t = /u 2 ) (57) 

with an arbitrary renormalization scale \x. At tree level, A(s, t) is momentum independent and with the 
proper normalization we have indeed X r (/J,) = A. 

Beyond tree level, the amplitude has an ultraviolet divergence that we regularize with a simple 
momentum cutoff A here. The relevant diagrams up to one loop are shown in Fig. 8. 

n 






+ ^ K + 



Fig. 8: Scattering amplitude for 00 — * <fi<j) to one-loop order. 

Setting s = — t = /j, 2 , one finds ((3q is a constant) 

X r (p) = A(s = -t = n 2 ) = A + /? A 2 logA/^ (58) 

+ [/,— independent terms of 0{\ ) + 0(\ ) . 

Since A r (/x) is finite, being equal to the physical scattering amplitude at some fixed point in phase space, 
the bare coupling A diverges as the cutoff A — ► oo. However, the bare coupling is not related to any 
physical quantity. Therefore, we are free to "sweep the infinities under the rug" as long as this is done in 
a transparent and controllable way. 

To do this, we change the renormalization scale by a small amount Sfi: 

X r ( P + Sn) - X r (p) = /3 A 2 log (—^-£\ + 0(A 3 ) 

= /3 A 2 log — ^— + 0(A r 3 ) = -/3 A 2 ^ + 0[(^) 2 ] + 0(A 3 ) . (59) 

The bare coupling and the cutoff have disappeared in the last equation. Expanding A r (/z + S/S) around fi 
and letting Sfi —> 0, we recover the (3 function of <p 4 theory to one-loop order: 

^dKiri = _ /3qA 2 (/x) + 0(A 3 } = p {Xr (n)) . (60) 

Unlike in Yang-Mills theories, (3q < so that 4> A theory is ultraviolet unstable like all quantum field theo- 
ries except nonabelian gauge theories [5, 12]. However, for understanding the essence of renormalization 
the important observation is that physical quantities do not depend on the bare coupling constant A nor 
on the cutoff A (for A —>■ oo) but only on the renormalized coupling A r (//) . For the purpose of comparing 
theory with experiment at present energies, we will never notice the stuff that was swept under the rug. 

We now turn to the choice of a regularization scheme. Although there are infinitely many possi- 
bilities, some choices are clearly better than others. The main criteria are: 

'•" The regularization method should respect symmetries of the theory as much as possible. In this re- 
spect, the previously employed momentum cutoff is as bad as it gets violating Poincare symmetry, 
gauge invariance, etc. 

'•" The scheme should violate only those symmetries that are necessarily violated by quantum effects 
(anomalies). 

*" The method should be simple to handle in practice. 

From the practitioner's point of view, dimensional regularization is the almost unique choice fulfilling 
these criteria. Let us demonstrate the method with a simple example, electronic vacuum polarization 
(setting m e = 0). 
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Gauge invariance, guaranteed by dim. regularization, 
implies the same structure as in the hadronic case: 

k-q n^(g) = (-f« 2 + ^)ntf) (61) 

8e 2 r(e) f 1 dxx(l-x) 



Fig. 9: Vacuum polarization at one loop. ^ ' (4-7r) 2 ~ e / n I— g 2 x(l — x)] e 

T(x) = l/x-7 + 0(x), 2e = 4-d. 

Since dimensional regularization works in d dimensions, there is a small problem here: there is no scale 
for log (— q 2 ) that will appear in the explicit form of H(q 2 ). To solve the problem, we insert unity (the 
scale \i and the constant c are completely arbitrary) in the expression and expand the second factor in e: 

l = (cfi)- 2e (cfi) 2£ = (cfi)- 2e [l + e log fi 2 + 2e log c + 0(e 2 )] . (62) 

Various schemes on the market differ by the constant c: 

MS c= 1 

MS logc = (7-log47r)/2. 

Using the most popular scheme (MS), the final result is 

I% 2 ) = -^ (M^ _ log (-^ 2 ) + |U 0(e) 



= uT(e^) ~ yj^ {log(-g 2 /^ 2 ) " §} • (63) 

The divergent part 11^ (e, jj) has been isolated and it will be absorbed by wave function renormalization 
of the photon field contributing to charge renormalization. We also notice that the coefficients of 1/e and 
— log (—q 2 //j, 2 ) are identical: the (3 function can be extracted from the divergent part. The stuff under 
the rug is useful after all. 

Back to the loop corrections of Fig. 7, we observe that the diagrams give rise to an amplitude 
proportional to g 2 , whereas A(e + e~ — ► qqG) ~ g s . How can the infrared divergences cancel among 
amplitudes of different order in g s ? The answer is that they cannot cancel on the level of amplitudes 
because the final states are different. Not the amplitudes but the rates must be added. Interference with 
the tree amplitude produces an 0(a s ) term in a(e + e~ — ► qq) that can and will cancel the infrared 
divergence in a(e + e~ — ► qqG). For details of the calculation I refer to the monograph [13], an excellent 
source for applications of perturbative QCD in general. 

The easier part are the loop corrections for a(e + e~ — ► qq). With dimensional regularization to 
regularize the infrared divergences, one finds 

^interference = ^^gtf^ /_! _ £ _ 16 + ( £ ) l > H{Q) = X . (64) 

The less familiar part is the three-body phase space integration in d dimensions giving rise to 

a ( 4 6 1 

(Tq qG = a C F ^H(e) I ^ + - + 19 + 0(e) | . (65) 

The infrared divergences cancel as expected. 



Adding the lowest-order cross section, one obtains finally 

cr(e + e~ -► hadrons) = a (l + 3CV j 1 + 0(a 2 )) = a (l + — + 0(a 2 )) , (66) 

with ctq defined in Eq. (54). The cross section cr(e + e~ — ► hadrons) is nowadays known up to 0{cfi s ). 
Replacing a s by the renormalization group improved running coupling a s (\/s), the general result for 



R e + e - can be written 



R P 



N c Y,Q)/e 2 
f 



1 + ^ C (^E>) 



n>l 



R 



(0) 



1 + Ci^U 



7T 



7T J 



A) 



C 2 -C^ log (s/f, 2 ) 



a s (/j,) 



TX 



+ 



(67) 



The normalization is such that C\ = 1, the coefficients C*2, C3 being also known. 

In principle, R e + e - (s) is independent of the arbitrary scale \x by construction. In reality, the 
unavoidable truncation of the perturbative series introduces a scale dependence. Although there is no 
unique prescription for the optimal choice of fj,, the obvious choice here is fi 2 = s to avoid large log- 
arithms. If the perturbative expansion is to make sense, we expect higher orders to mitigate the scale 
dependence. This is nicely demonstrated in Fig. 10 taken from Ref. [13] where the deviation (in percent) 
of R e + e - (y/s = 33 GeV) from R+ _ is plotted as a function of fi. As expected, ft 2 = s is indeed a 
very reasonable choice already at 0(a s ) (L). 
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Fig. 10: Improvement of the scale dependence in higher orders of perturbation theory for R e + e - (\/s = 33 GeV) 
(taken from the book of Ellis, Stirling and Webber [13]). 



2.3 Measurements of a s 

How should one characterize the coupling strength of QCD? After all, the scale \x is arbitrary and, in 
addition, a s {jx) is in general scheme dependent. For qualitative purposes, one may introduce a scale 
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Aqcd that is independent of the renormalization scale fj,. The drawback is that this quantity is scheme 
independent only at leading (one-loop) order: 

4tt 



a s {E) 



/?olog(£7A^ C D)' 
The coefficient j3q is defined by rewriting the [3 function for a s (instead for g s as in Sec. 1.1): 






Po_ 2 _ Jh_ 3 , 

2vr s 4tt 2 s + '" 



(68) 



(69) 



The j3 function is known up to four loops (coefficient j3z) but only the first two coefficients 

fa = 11 - 2N F /3 , fa = 51 - 19N F /3 (70) 

are scheme and gauge independent. 

Since the scheme dependence is unavoidable, the coupling strength is nowadays usually given in 
the form of af s (Mx). Of course, this fixes a^ s (fi) at any scale via the integral 

ra s (M2) dx 



2 /..2\ 



log (/^2/V 



,(W) P( x ) 



To get a first rough estimate of a s , consider the leading-order prediction 



R e+e -(Mz) = RfL(M z )[l + 



a s (M z ) 



7T 



(71) 



(72) 



Comparing the combined LEP result [14] R e + e - (Mz) = 20.767 ± 0.025 with the tree-level prediction 
R^ e - {M z ) = 19.984, one obtains 

a s {M z ) = 0.123 ±0.004, (73) 

close to the full three-loop result and not bad at all for a first estimate. A compilation of results can be 
found in the Review of Particle Properties [14]. In Figs. 11 and 12, the most recent data compiled by 
Bethke [15] are shown. 
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Fig. 11: Compilation of data for the extraction of 
a^iMz) by Bethke [15]. 
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Fig. 12: Energy dependence of the running coupling con- 
stant ap(Q) [15]. 
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The first impression is the remarkable agreement among experiments and with theory. However, for 
determining the best value of a s (Mz), the following two problems must be kept in mind: 

- Different observables are known with different theoretical accuracy: next-to-leading order (NLO) 
vs. next-to-next-to-leading order (NNLO). Not only the scale dependence but also different scheme 
dependences must be taken into account. 

- Theoretical errors are not normally distributed. 

Using only NNLO results, Bethke found [15] 

a s (M z ) = 0.1182 ± 0.0027 , (74) 

very similar to the PDG average [14] (using a different procedure) 

a s (M z ) = 0.1187 ±0.0020. (75) 

All values in this paragraph refer to a^ s (Mz)- 

2.4 Hadronic r decays 

A remarkably precise value for a^ s (Mz) comes from hadronic r decays. At first sight, this is quite 
surprising because at the natural scale fi = m T one has approximately a s {m T ) ~ 0.35. Can one ex- 
pect reasonable convergence of the perturbative series for such a large coupling and how big are the 
nonperturbative corrections? 

The first systematic investigation of R T = r(r~ —> u T + hadrons)/r(r~ — ► v T e~V^,) was 
performed by Braaten, Narison and Pich [16]. The analysis is similar to the one for R e + e - , with obvious 
modifications: the electromagnetic current (coupling to e + e~) must be replaced by the charged weak 
current (coupling to tu t ). 

We start again with the two-point function (of weak currents L M = «7 M (1 — 75)^0): 
nf(g) = i j d A xe i( i' x {Q\TL^{x)L v {^\Q) 

= {-gTf + gV) n? } (q 2 ) + gV4 0) (q 2 ) , (76) 

with dg the Cabibbo rotated d-quark field. One major difference to the electromagnetic case is that one 
has to integrate over the neutrino energy or, equivalently, over the hadronic invariant mass s: 

fl ' =i2 'r$( i -^) 2 {( i+2 ^) imn ™ (s)+imn ™ <s) } • <77> 

The problem is that the integration extends all the way down to s = (for m u = nid = 0) where pertur- 
bation theory is certainly not applicable. 

However, QFT provides information about the analytic structure of two-point functions that can 
be used in a standard manner to circumvent the problem. The invariant functions n^ ' (s) are known to 
be analytic in the complex s-plane with a cut on the positive real axis. Therefore, Cauchy's theorem tells 
us that the contour integral in Fig. 13 vanishes. 

One can now trade the integral along the cut of 



r(0,l), 



1 



rC ' 1 )/,, ■ .vi-ttC . 1 ), 



ImU^(s) = - [Wis + is) - nr> " is) 



(78) 

Vi L _ - J 

for an integral along the circle \s\ = m^ in the complex s-plane. It turns out that the nonperturbative 
corrections are now manageable, being suppressed as (Aqcd/^v) 6 - Very helpful in this respect is the 

( s \ 2 
factor 1 7T in the integrand that suppresses potentially big contributions near the endpoint of the 

cut. 
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Fig. 13: Contour in the complex s-plane for the two-point functions IT^ ' (s). 
The final result can be written in the form [16] 

Rr = 3(\Vud\ +\V US \ ) 'S'EW { 1 + ^EW + ^pcrt + ^nonpcrt } (79) 

with leading and nonleading electroweak corrections Sew = 1-0194 and c^ w = 0.0010, respectively. 
The perturbative QCD corrections of interest for the extraction of a s are contained in 

a s (m T ) f 19 \ ( a s (m T ) \ 2 

<W - —^ + \C 2 + -fioj [^^) +■■■ (80) 



IT \ IT J \ TT 

Finally, the best estimates of nonperturbative contributions, using QCD sum rules and experimental input, 
yield 

^nonpcrt = "0.014 ± 0.005 . (81) 

From the PDG fit for R T one then obtains a s {m T ) = 0.35 ± 0.03. More interestingly, running this value 
down to Mz with the help of the four-loop (5 function, one finds 

a s (M z ) = 0.121 ± 0.0007(exp) ± 0.003(th), (82) 

not only compatible but in fact very much competitive with other high-precision determinations. 

2.5 Deep inelastic scattering 

From the conception of QCD till today, deep inelastic scattering of leptons on hadrons has had an enor- 
mous impact on the field. It is also a classic example for the factorization between long- and short- 
distance contributions. 

Let us start with the kinematics of (in)elastic electron-proton scattering e~(k) + p(p) — ► e~ (k r ) + 
X(px) shown in Fig. 14. In the case of elastic scattering (X = p), we have 

W 2 = m 2 , Q 2 = 2mv , x = l (83) 

and the usual two variables are s = (p + k) 2 , Q 2 with the differential cross section da(s, Q 2 )/dQ 2 . 

For the inclusive scattering, there is a third independent variable: a convenient choice is s, x, y 
with < x < 1, < y < 1. In general, we distinguish different types of deep inelastic scattering: 

1Q 



type exchange 

neutral current (NC) DIS 7, Z, 7Z-interference 

charged current (CC) DIS W ± 

For these lectures, I will restrict the discussion to photon exchange and to unpolarized (spin-averaged) 
DIS. 
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W 2 =p 2 x = ip + q) 2 = m 2 + Imv - Q 2 > m 2 . 



Fig. 14: Deep inelastic scattering. 
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The matrix element for the diagram in Fig. 14 has the structure eZ(epton)^ —r e /i(adron);,, with 

leptonic and hadronic current matrix elements /(epton) M and ^(adron)^, respectively. The resulting 
double differential cross section is of the form 

d 2 <r , 9, d 2 a 2-Kiia 2 „ TT „„ 

i- 1 ^ = 2[k^k v ' k^kv — k ■ k g^y) 
H^(p, g) = i- | d*ze*1 ■ z (p, s\ [j»Jz], J e V(0)] \p, s) . 

One now performs a Lorentz decomposition of the hadronic tensor H^ v and contracts it with the leptonic 
tensor L^ (setting m e = 0). The differential cross section then depends on two invariant structure 
functions i<i, i*2 (in the photon case), which are themselves functions of the scalars p.q, q 2 or v, Q 2 or 
x,Q 2 : 

d 2 o Q 2 d 2 a Wf^^W ^ ^ ^1 (g5) 



dxdy y dxdQ 2 xyQ 2 \\ " Q 

Deep inelastic scattering corresponds to Q 2 3> m 2 and W 2 S> m 2 . While the cross section shows a 
rather complicated behaviour at low and intermediate momentum transfer, the structure functions exhibit 
an originally unexpected simple behaviour in the so-called Bjorken limit 

Q S> m , v S> m with x = fixed . (86) 

2mv 

As shown for F2 in Fig. 15, in the Bjorken limit the structure functions seem to depend on the variable x 
only: 

Fi(x,Q 2 ) -^ Fi(x) . (87) 

This scaling behaviour suggested that the photon scatters off point-like constituents (no scale) giving rise 
to the quark parton model (QPM). 



?() 




Fig. 15: Evidence for Bjorken scaling taken from Ref. [13]. 



QPM in the Breit frame 



The characteristics of the QPM can best be visualized in the so-called Breit frame where the proton and 
the virtual photon collide head-on. 



Breit frame 



vAAAAAAAA^ 

q 



p ~ (P, P, 0, 0) with P > m 
q = (0,-^Q2,0,0) 



Fig. 16: DIS in the Breit frame. 

The nucleon is pictured as a bunch of partons with negligible transverse momenta. Each parton carries a 
fractional momentum £p. Since the scattered quark is massless (compared to P), we have 



(q + ivf 



and therefore 



£ 



P 



V 2 



2x 



-Q 2 + 2&.q = 



q + xp=(xP,-VQ^/2,0,0) . 



(88) 



(89) 



The struck parton scatters with momentum q + xp backwards, i.e. in the direction of the virtual photon, 
justifying a major assumption of the QPM: the virtual photon scatters incoherently on the partons. 

The fundamental process of the QPM is elastic electron-quark scattering 



e (k) + q(£p) — ► e (k') + q(£p + q) ■ 



(90) 



Since there are now only two independent variables, the double differential cross section in Eq. (85) 
contains a 5-function setting x equal to £: 



d 2 a 



(?) 



dxdy yQ 



^[i + d-y)T^(*-0. 



(91) 
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In the notation of Eq. (85), 



F 2{q) = xQ 2 q 5(x - = 2xF 1 



(9) • 



(92) 



The incoherent sum of partonic cross sections amounts to an integral over quark distribution functions 



Q,Q ° 1,9 



implying the Callan-Gross relation [17] 



F 2 (x) =2xF 1 (x) 



(93) 



(94) 



that is due to the spin- 1/2 nature of quarks. The so-called longitudinal structure function Fl = F 2 — 2xFi 
therefore vanishes in the QPM. 

It was already known at the beginning of the seventies, before the advent of QCD, that exact 
scaling in the sense of the QPM was incompatible with a nontrivial QFT. QCD must therefore account 
for the systematic deviation from scaling that is clearly seen in the data (e.g., in Fig. 17): with increasing 
Q 2 , the structure function F 2 increases (decreases) at small (large) x. 
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Fig. 17: Experimental evidence for the violation of Bjorken scaling taken from Ref. [18]. 
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Qualitatively, scaling violation is due to the radiation of (hard) gluons generating transverse momenta for 
the quarks. More gluons are radiated off when Q 2 increases, leading to logarithmic scaling violations in 
the structure functions and to scale dependent partem distribution functions (pdf) qi(x,fi 2 ), g(x,fi 2 ) as 
we will now discuss in more detail. 



DIS in QCD 



In QCD at leading order in g s , the same diagrams need to be considered as the ones in Figs. 6, 7 relevant 
for e + e~ — ► hadrons except for crossing. 

Previously, the sum of real and virtual gluon emission was infrared finite because a sum over all 
final state quarks and gluons was involved. In DIS the situation is different because the initial state con- 
tains a quark. Since different incoming quark momenta are in principle distinguishable, gluons collinear 
with the incoming quark generate in fact an infrared divergence. How to get rid of this divergence will 
be discussed later but for now we regulate the infrared divergence with a cutoff on the transverse quark 
momentum k\ > k 2 . Adding the contribution to F 2 ( q ) from real gluon emission (£ = 1 for simplicity) 
one finds 

a, („ , ,, Q 2 



F 2{q) (x,Q z ) = Q z q x 



5(1 



2tt 



+ -±(P qq ( x )log^ + C q ( 



(95) 



Introduction of the cutoff k has produced a logarithmic dependence on Q 2 . This dependence is governed 
by the so-called quark-quark splitting function P qq (x). This function is universal, in contrast to the 
non-logarithmic coefficient function C q (x) that is scheme dependent. 

The origin of the scheme-independent log Q 2 can be understood as follows. The struck quark 

dk 2 
acquires a transverse momentum k± with probability a s — ^. Since k 2 , cannot be bigger than Q 2 , 

kj_ 

integrating over all k± produces the term a s log Q 2 /k 2 . Virtual gluon contributions (diagrams in Fig. 7) 
must still be added and they are ultraviolet finite as before. Since this is a contribution from elastic 
scattering it must be proportional to 6(1 — x). Altogether, the quark-quark splitting function at leading 
order is 

p -w=Kra) +2i(1 ^'' <96) 

Also the first part of P qq (x) is actually a distribution. The distribution [F(x)] + is defined in such a way 
that for every sufficiently regular (test) function f(x) one has 



dxf(x)[F(x)l 



dx(f(x)-f(l))F(x) 



(97) 



It is then straightforward to show that P qq (x) in Eq. (96) can also be written as 



P qq (x) 



1 + x 2 
3 [(1 — a?) 



(98) 



The problem remains how to interpret (or rather get rid of) the infrared cutoff k. Up to now, we 
have only considered the quark structure function F 2 ( q )(x, Q 2 ). To get F 2 (x, Q 2 ) for the nucleon, we 

convolute F 2 ^(—,Q 2 ) with a (bare) pdf go(0 : 



F 2 (x, Q 2 ) = x^2 Q\ \qo(x) + J f %o(y) [P qq (x/v) log J + C q (x/y) 

q>q l n Jx y i 



(99) 
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One now absorbs the collinear singularity 
renormalized pdf q(x, ji 2 ): 



logK 2 into qo(x) at a factorization scale /j, to define a 



q(x,fi 2 ) 



qo(x) + g J 1 % (y) {P qq (x/y) log ^ + C' q (x/y)\ . (100) 



The interpretation of q(x, fi 2 ) is straightforward: the soft part k\ < [i 2 is now included in the pdf. Since 
the scale \x is arbitrary, the renormalized pdf is necessarily scale dependent, in complete analogy with the 
renormalization of ultraviolet divergences. As a small aside, we take note that the coefficient function 
C' q need not be the same as C q in Eq. (95), both being scheme dependent. The final form for the nucleon 
structure function F 2 in the MS scheme (except for a contribution from the gluon pdf) is then 
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This factorization formula can be proven to all orders in a s , separating the calculable hard part from 
the soft part contained in the scale dependent pdfs. The pdfs q(x,fi 2 ), q(x,^ 2 ), g(x,^ 2 ) describe the 
composition of nucleons and are, of course, not calculable in perturbation theory. They can be extracted 
from experimental data with the help of appropriate parametrizations (cf., e.g., Ref. [13]) but the question 
remains what the factorization result (101) is actually good for. 

The answer is that even though the functional dependence of the pdfs cannot be calculated their 
scale dependence is calculable in QCD perturbation theory. The derivation of the so-called DGLAP 
evolution equations [19] is very similar to the derivation of the (3 function in Sec. 1.1, starting from the 
observation that the measurable structure function must be scale independent: 
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With Pq q (x) given by Eq. (98), the evolution equation at leading order takes the explicit form 



H 



<d q(x,iJ, 2 

d[i 2 



2« s (^) 
3vr 



f 1 dz 2 1 

/ —q{x/z,fi )- 

Jx Z 1 



1 + z 2 



2a s (^) 
3vr 



q(x,n ) I dz 

o i 



1+z 2 



(104) 



Due to soft gluons, both terms on the right-hand side are divergent: the first term with positive sign is 
due to quarks with momentum fraction larger than x radiating off gluons whereas the second term leads 
to a decrease from quarks with given x that radiate gluons. The overall result is finite. 

At 0(a s ) also the gluon pdf enters via 7* + g —> q + q. At any order, the DGLAP equations are 
in general (2^^ + 1) -dimensional matrix equations for qi(x,fi 2 ), 7ji(x,[i 2 ) (i = 1, . . . ,Nf), g(x,fi 2 ) 
with splitting functions P ( 



qq 



[x), Pqg(x), Pgq(x) and P gg (x). The analytic calculation of these splitting 
functions to next-to-next-to-leading order (three loops) has just been completed [20] allowing for precise 
tests of scaling violations. Finally, we note that the longitudinal structure function Fl(x, Q 2 ) that was 
zero in the QPM is generated in QCD already at 0(a s ). 

For more applications of perturbative QCD I refer once again to the book of Ellis, Stirling and 
Webber [13]: jets in e + e _ and hadroproduction, vector boson production (Drell-Yan), heavy quark pro- 
duction and decays, Higgs production at the LHC, . . . 
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3 Heavy and light quarks 

3.1 Effective field theories 

Unlike QED, QCD is valid down to shortest distances because of asymptotic freedom. However, at long 
distances where quarks and gluons are practically invisible, perturbative QCD is not applicable and a 
nonperturbative approach is needed. Many models can be found in the literature that are more or less 
inspired by QCD. Qualitative insights into the structure of the strong interactions have been found from 
model studies but quantitative predictions require methods that can be related directly to QCD. There are 
essentially only two approaches that satisfy this criterion. 

w Lattice QCD has already scored impressive results and may in the long run be the most predictive 
method. At present, the range of applicability is still limited. 

w Effective field theories (EFTs) are the quantum field theoretical formulation of the "quantum lad- 
der": the relevant degrees of freedom depend on the typical energy of the problem. EFTs become 
practical tools for phenomenology when the characteristic energy scales are well separated. 

Let a given step of the quantum ladder be characterized by an energy scale A. The region E > A is the 
short-distance region where the fundamental theory is applicable. At long distances (E < A), on the 
other hand, an effective QFT can and sometimes must be used. By definition, the notions "fundamental" 
and "effective" only make sense for a given energy scale A. As we probe deeper into the physics at short 
distances, today's fundamental theory will become an effective description of an even more fundamental 
underlying theory. 

To understand the different effective field theories that are being used in particle physics, it is 
useful to classify them as to the structure of the transition between the fundamental and the effective 
level. 

• Complete decoupling 

The heavy degrees of freedom (heavy with respect to A) are integrated out, i.e., they disappear 
from the spectrum of states that can be produced with energies < A. Correspondingly, the effective 
Lagrangian contains only light fields (once again, light stands for masses < A) and may be written 
symbolically as 

£eff = £d<4 + Yj JdTZ E 9i d O id . (105) 

d>4 i d 

The first part Cd<i contains all the renormalizable couplings for the given set of fields. The 
best-known example for such a Lagrangian is the Standard Model itself. The second part of the 
Lagrangian (105) contains the nonrenormalizable couplings having operator dimension d > 4. 
The best-known example here is the Fermi theory of weak interactions with d = 6 and A = M\y. 
For the Standard Model, on the other hand, we do not know the scale where new physics will 
appear. Present experimental evidence implies A > 100 GeV but there are good reasons to expect 
new physics around A ~ 1 TeV. 

• Partial decoupling 

In this case, heavy fields do not disappear completely in the EFT. Via so-called field redefini- 
tions, only the high-momentum modes are integrated out. The main application of this scenario in 
particle physics is for heavy quark physics. 

• Spontaneous symmetry breaking (SSB) 

In the previous two classes, the transition from the fundamental to the effective level was smooth. 
Some of the fields or at least their high-energy modes just drop out and the effective description 
involves the remaining fields only. In the present case, the transition is more dramatic and involves 



a phase transition: SSB generates new degrees of freedom, the (pseudo-)Goldstone bosons associ- 
ated with spontaneously broken symmetries (to be discussed in more detail in Sec. 3.4). The prefix 
pseudo accounts for the frequent case that the symmetry in question is not only spontaneously but 
also explicitly broken. Goldstone bosons in the strict sense are massless and the associated SSB 
relates processes with different numbers of Goldstone bosons. As a consequence, the distinction 
in the Lagrangian (105) between renormalizable (d < 4) and nonrenormalizable (d > 4) terms 
becomes meaningless. Therefore, EFTs in this category are generically nonrenormalizable. An 
important but maybe too simple exception is the Higgs model for electroweak SSB. 

The generic EFT Lagrangian is organized in the number of derivatives of Goldstone fields and in 
the number of terms with explicit symmetry breaking. An important concept is universality: it 
turns out that EFTs describing different physical situations have very similar structure. In QCD, 
the symmetry in question is chiral symmetry that becomes exact in the limit of massless quarks. 
In the real world, SSB of chiral symmetry generates pseudo-Goldstone bosons that are identified 
with the pseudoscalar mesons it, K,rj. 

We are used to derive quantitative predictions from renormalizable QFTs in the framework of perturba- 
tion theory but how should we treat nonrenormalizable EFTs? The clue to the answer is Goldstone's 
theorem [21] that makes two crucial predictions. The first prediction is well known: SSB implies the ex- 
istence of massless Goldstone bosons. The second consequence of Goldstone's theorem is not that well 
known but very important as well: Goldstone bosons decouple when their energies tend to zero. In other 
words, independently of the strength of the underlying interaction (the strong interaction in our case!), 
Goldstone bosons interact only weakly at low energies. This important feature allows for a systematic 
expansion of strong amplitudes even in the confinement regime, which is precisely the low-energy regime 
of QCD. 

However, in contrast to the decoupling case (e.g., in heavy quark physics), the coupling constants 
of the low-energy EFT cannot be obtained by perturbative matching with the underlying theory of QCD. 
Other methods have to be used to get access to the low-energy couplings. 

3.2 Heavy quarks 

Quarks cannot be put on a balance, or more realistically, their energies and momenta cannot be measured 
directly. How do we then determine their masses? Two methods have been used. 

<•" The first approach ignores confinement and calculates the pole of the quark propagator just as 
we determine at least in theory the mass of the electron. This looks rather artificial because the 
full quark propagator should have no pole because of confinement. Going ahead nevertheless, one 
expects those pole masses to be very much affected by nonperturbative infrared effects. In practice, 
this method is only used for the top quark with [14] 

m t = 174.3 ±5.1 GeV. (106) 

*" Quark masses are parameters of the QCD Lagrangian j ust like the strong coupling constant g s . One 
therefore studies the influence of these parameters on measurable quantities and extracts specific 
values for the masses by comparison with experimental measurements. As for the strong coupling 
constant, renormalization is crucial and introduces a scale dependence also for quark masses. The 
scale dependence is governed by a differential equation very similar to the renormalization group 
equations (9) or (69) for the strong coupling: 

^dm^l = _ tMm)) (//) (107) 
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where the anomalous dimension 7 is nowadays known up to four-loop accuracy: 
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The solution of this renormalization group equation for m q (//) is flavour independent (in the MS 
scheme): 
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Since 7(a s ) is positive, quark masses decrease with increasing scale \x, e.g., 

m,(l GeV) 



(109) 



mJM. 



2.30 ±0.05. 



(110) 
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Different methods have been used to determine the quark masses: H(eavy) Q(uark) E(ffective) 
T(heory) (see below), QCD sum rules (section 3.3), lattice QCD, . . . The current state is summarized in 
Fig. 18 taken from Ref. [22]. All values correspond to the MS scheme: light quarks are given at the scale 
2 GeV whereas the heavy quarks m c , nib are listed as m q (m q ). 
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Fig. 18: The values of quark masses taken from the 2004 Review of Particle Properties [22] . The most recent data 
points are at the top of each plot. 
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Heavy quark effective theory (HQET) 



Why should one use an EFT for 6-quark physics? After all, QCD is still accessible in perturbation theory 
for \x = m&. The main arguments in favour of HQET are the following. 

•" Although the hard effects are calculable in QCD perturbation theory, there are inevitably incal- 
culable soft effects because hadrons rather than quarks and gluons appear in the final states of B 
decays. The necessary separation between perturbative and nonperturbative contributions is much 
easier to achieve in an EFT description. The keyword is the same as in deep inelastic scattering: 
factorization. 

•* Approximate symmetries that are hidden in full QCD become manifest in an expansion in 1/mg. 

••" Explicit calculations simplify in general, in particular the resummation of large logs via renormal- 
ization group equations. 

The physics behind HQET can be understood by an analogy with atomic physics. The spectrum of the 
hydrogen atom is to a good approximation insensitive to the proton mass. In fact, the atomic spectra of 
hydrogen and deuterium are practically identical. The implementation of this analogy is most straight- 
forward for hadrons with a single heavy quark (b or c). In the hadron rest frame the heavy quark "just 
sits there" acting as a colour source just as the proton acts as a Coulomb source in the hydrogen atom. 

We decompose the momentum of the heavy quark as 

p^ = m Q v^ + k^, (111) 

where v is the hadron velocity normalized to v 2 = 1 (v = (1, 0, 0, 0) in the hadron rest frame). The 
residual quark momentum k is then expected to be of O(Aqcd) only. Starting from the QCD Lagrangian 
for a heavy quark Q, 

C Q = Q(0 - m Q )Q , (112) 

we decompose the quark field Q(x) into two fields h v (x) and H v (x) by using energy projectors P^ 1 = 
(1 ± f)/2 and applying (shift) factors e tm Q v - x -_ 

Q( x ) = e -im Q vx (/ ^ x) + Hv ( x ^ (U3) 

h v {x) = e im ^ v - x P+Q(x) , H v {x) = e im ^ x p-Q(x) . 

It is easy to check that in the hadron rest frame the fields h v and H v are just the upper (big) and lower 
(small) components of the spinor field Q, respectively. Expressing Cq in terms of h v and H v , one finds 

C-Q = Q{i$ - ™q)Q 

= h v iv ■ D h v — H v (iv ■ D + 2m,Q)H v + mixed terms . (114) 

For the purpose of illustration, we use the field equation (ilj) — rciQ)Q = to eliminate H v in this 
Lagrangian: 

C Q = KivDh v + Kilp 1 _ 0±h v with D>\ = (g^ - v»v v )D v . (115) 

iv ■ D + 2toq — ie 

The heavy quark mass itlq has disappeared from the kinetic term of the shifted field h v and has moved 
to the denominator of a nonlocal Lagrangian that is in fact the starting point for a systematic expansion 
in 1/rriQ. 

The Lagrangian for b and c quarks to leading-order in 1/rriQ is therefore 

A>,c = b v iv ■ Db v + c^iv ■ D c v . (116) 



This Lagrangian exhibits two important symmetries. The symmetries are only approximate because the 
Lagrangian (116) is not full QCD but the first approximation in an expansion in 1/tbq. The symmetries 
are manifest in (1 16) but they are hidden in full QCD. 

w A heavy-flavour symmetry SU(2) relates b and c quarks moving with the same velocity. 

w Because there is no Dirac matrix in the Lagrangian (1 16), both spin degrees of freedom couple to 
gluons in the same way. Together with the flavour symmetry, this leads to an overall spin-flavour 
symmetry SU(4). 

The simplest spin-symmetry multiplet M. consists of a pseudoscalar M and a vector meson M*. One of 
the first important applications of spin-flavour symmetry was for the semileptonic decays B — ► D^lui. 
In general, there are several form factors governing the two matrix elements (for D and D*). To leading 
order in 1/rriQ, all those form factors are given up to Clebsch-Gordan coefficients by a single function 
£(v • v') called Isgur-Wise function: 

(M(v')\h^Th v \M(v)) ~ £(v ■ v') . (117) 

Moreover, h v ^h v is the conserved current of heavy-flavour symmetry. Similar to electromagnetic form 
factors that are normalized at q 2 = due to charge conservation, the Isgur-Wise function is fixed in the 
no-recoil limit v = v' to be 

£(vV = l) = l. (118) 

Of course, there are corrections to this result valid only in the symmetry limit, both of 0(a s ) and in 
general of 0(l/mg). For the decay B — ► D*lv\, the leading mass corrections turn out to be of 0(1/171%) 
only. HQET provides the standard method for the determination of the CKM matrix element Vd, (see 
Ref. [23] for reviews). 

HQET has been extended in several directions. 

• Soft collinear effective theory (SCET) 

HQET cannot be applied to decays like B — ► X s j or B — ► -kit where the light particles in the 
final state can have momenta of O(tbq). SCET accounts for those energetic light states but it is 
more complicated than HQET. Because of the presence of several scales, several effective fields 
must be introduced by successive field transformations. A major achievement is again the proof 
of factorization that is for instance crucial for a reliable extraction of the CKM matrix element V u j, 
from experiment. 

• Nonrelativistic QCD (NRQCD) 

This extension of HQET includes quartic interactions to treat heavy quarkonia b b and c c. In this 
case, three widely separate scales are involved: the heavy mass uiq, the bound-state momentum 
p ~ itiqv (v <C 1) and the kinetic energy E ~ itiqv 2 . Applications include the analysis of the q q 
potential and the production and decay of quarkonia [24]. 

3.3 QCD sum rules 

The general idea of QCD sum rules is to use the analyticity properties of current correlation functions 
to relate low-energy hadronic quantities to calculable QCD contributions at high energies. We recall the 
example of the two-point functions n^ ' (q 2 ) in hadronic r decays discussed in Sec. 2.4. 
In general, the QCD contribution consists itself of two different parts, 

- a purely perturbative part and 

- a partly nonperturbative part that is important at intermediate energies and makes use of the oper- 
ator product expansion (OPE, another case of factorization). 



Altogether, a typical two-point function (QCD sum rules are not restricted only to two-point functions, 
however) has the form 

U(q 2 ) =Il pcrt (q 2 ) +Y,C d (q 2 )(0\O d \0) . (119) 

d 

The so-called Wilson coefficients Cd(q 2 ) are calculable perturbatively and they depend on the two-point 
function under consideration. Up to logs, they decrease for large \q 2 \ as (g 2 ) _n w with some positive 
integer n(d). The vacuum condensates (0|O<f |0), on the other hand, are universal and they absorb long- 
distance contributions with characteristic momenta < \/jg^|. 

Three main types of applications of QCD sum rules can be distinguished. 

• Using experimental data as input for the low-energy hadronic part, one can extract universal QCD 
parameters: a s , quark masses, condensates, . . . 

• With QCD parameters known, one can predict hadronic quantities: hadron masses, decay con- 
stants, amplitudes, . . . 

• The compatibility of low-energy data with QCD can be checked. I will discuss a recent example 
of topical interest, the spectral data relevant for the leading hadronic contribution to the anomalous 
magnetic moment of the muon. 
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The biggest source of uncertainty in the calculation of the anomalous magnetic moment of the muon 
<V = (ST* ~~ 2)/2 in the Standard Model is at present the lowest-order hadronic vacuum polarization 

/>oo 

had,LO = a vac. P oi. = / dtK(t)a (e + e~ -» hadrons)(t) (120) 

depicted in Fig. 19. 
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Fig. 19: Lowest-order hadronic vacuum polarization contribution to the anomalous magnetic moment of the muon. 

The kernel K(t) is a known function [25]. Although the integral extends from threshold to infinity, 
about 73 % of a M a ' are due to the -k-k intermediate state in Fig. 19, governed by the pion form factor. 
Moreover, 70 % of a^ ' come from the region t < 0.8 GeV 2 . Therefore, by far the most important 
part is not calculable in QCD perturbation theory. 

A few years ago, Alemany, Davier and Hocker [26] suggested to use not only data from e + e~ — » 
7r + 7r - to extract the pion form factor but also from the decay r _ — ► tt q it~v t . In the isospin limit, it is 
straightforward to derive the relation 

a {e + e- -► vr + vr-)(t) = h(t)—± j f Tj - (121) 

30 



with a known function h{t) where t is the two-pion invariant mass squared. 

At the level of precision required for comparison with experiment (better than 1 % for a M a ' ), 
isospin violating and electromagnetic corrections are mandatory [27,28]. The status until recently was 
summarized by Hocker at the High Energy Conference in Beijing [29]. 

w There was a significant discrepancy between the r and e + e~ data, mainly above the p resonance 
region, as shown in Fig. 20. 
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Fig. 20: Comparison between e + e~ and r data for the pion form factor from Ref. [29]; plotted is the difference of 
the form factors squared normalized to the r data. 



<•" Isospin violation cannot explain the difference. 

<•" The e + e~ data from the KLOE experiment confirm the previous trend of the CMD-2 data although 
the agreement among the e + e~ data is not impressive. 

The widely accepted recommendation at the Beijing Conference was to ignore the r data until the origin 
of the discrepancy is understood [29]. 

The situation has changed both on the theoretical and on the experimental side. The theoretical 
clarification is due to Maltman [30] who checked the consistency of experimental data with QCD by 
investigating sum rule constraints for the two spectral functions relevant for the e + e~ and the r case, 
respectively: 

Pem(s) =lmllei m (s) and /9y =1 (s) =lmn L ,ud( s ) ■ 

By using a contour integral in the complex s-plane as in Fig. 13, with rr? T replaced by an a priori arbitrary 
so. one derives a so-called F(inite)E(nergy)S(um)R(ule) 



f s ° 1 / 

/ w(s)p(s)ds = ——— <p w(s)H(s)ds 

JO 27T J\ s \ =so 



(122) 



II(s), p(s) refer to the two cases considered (electromagnetic currents in the e + e case, charged currents 
for t) and w(s) is an analytic function (actually a polynomial) that will be chosen conveniently. 



II 



The right-hand side of the FESR can be estimated in QCD as exemplified by Eq. (119). The purely 
perturbative part is known up to af, with estimates of the O(af) contribution available. For not too small 
so the purely perturbative part dominates the right-hand side depending only on a s (Mz)- The d = 2 part 
is completely negligible in the r case and it depends only on the mass of the strange quark in the elec- 
tromagnetic case. The d = 4 contributions involve the quark condensates (0\m q qq\0)(q = u, d, s) (very 
well known from chiral perturbation theory, cf. Sec. 3.5) and the gluon condensate (OlasG" Ga U \0) 
(less well but sufficiently known from charmonium sum rules). For d > 6, the relevant condensates are 
practically unknown. However, by using the pinching trick (w(so) = 0) appropriately, Maltman elimi- 
nates the d = 6 OPE contributions to II (s) altogether. The negligible effect of d > 8 contributions can 
be checked by varying so- 

Turning now to the left-hand side of the FESR (122), Maltman uses the most precise experimen- 
tal data for the spectral functions: ALEPH (compatible with CLEO and OPAL) for p^ l {s), CMD-2 
for p em (s). As a first test, he fits a s (Mz) (keeping the remaining OPE input fixed) from the exper- 
imentally determined left-hand side of the FESR. The outcome is, quite independently of the weight 
function w(s) that is always chosen positive and monotonically increasing for < s < sq, that the fitted 
value of a s (Mz) is systematically lower than the high-energy determination dominated by LEP in the 
electromagnetic case while there is perfect agreement in the r case. This is a first indication that the 
electromagnetic spectral density is too low. 

A second test, largely independent of the value of a s (Mz), compares the slopes of the OPE parts 
and spectral integrals with respect to sq. The results are shown in Figs. 21, 22. 
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Fig. 21: Slopes with respect to so in the e + e~ case for a 
specific weight function wq(s). The dashed lines denote 
the central values for the OPE input and the solid lines 
indicate the error bands. The spectral integrals are shown 
for several points, error bars included. 



Fig. 22: Same as in the previous figure for the t data. 
Both figures are taken from Ref. [30]. 



The situation is similar as before. While the slopes differ between data and QCD by > 2.5 a in the 
electromagnetic case there is perfect agreement in the r case. The previous conclusion is reinforced: the 
e + e~ spectral data are systematically too low whereas the r data are completely consistent with QCD, 
both in normalization and in the slopes. The QCD sum rule tests clearly favour the r over the e + e~ data 
for the pion form factor. 
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The most recent development is again an experimental one. Only two months before the School 
new data on e + e~ — ► -k + tt~ were released by the SND Collaboration from Novosibirsk. Their results 
disagree with both CMD-2 and KLOE and they go in the right direction towards reconciling the e + e _ 
with the r data. The discrepancies between the three data sets in e + e~ annihilation remain to be under- 
stood. Based on r data for the 2ir and Air channels in the hadronic vacuum polarization, the calculation 
of the anomalous magnetic moment of the muon in the Standard Model [31] compares well with the 
measured value [32] : 

< P - aJ M = (7-6 ± 8.9) • 1(T 10 . (123) 

There is at present no evidence for a discrepancy between the Standard Model and experiment. 

3.4 Chiral symmetry 

By construction, QCD is a gauge theory with gauge group SU(3) C . However, the QCD Lagrangian 

1 N F 

£qcd = --tr(GVGH + Y, Qf (0 ~ "Vic) Qf (124) 

/=i 
possesses additional symmetries. As in QED, the theory is parity invariant because of the absence of 75 
(vector couplings only). Moreover, the coupling constant g s and the quark masses mj are real so that 
QCD conserves also CP, ignoring the so-called strong CP problem here. 

Are there still additional symmetries in the QCD Lagrangian (124)? To answer this question, we 
first have a look at the quark kinetic term only (with Np = 6 flavours): 

6 6 

£kin = i ^2 q f Jp qf = i ^2 {QfL $ Q/l + QfR Ipqfn} , (125) 

with chiral components 

1 1 

QL = -(1 - 75)9, QR = ^0-+75)q- (126) 

Since the qjp and the qjR do not talk to each other in (125), they can be rotated separately implying the 
maximal global flavour symmetry U(6) x U(6). However, this is a symmetry of the kinetic term only. 
In the full quark Lagrangian (colour indices will be suppressed from now on) 

6 

£q = J2 {QfL *& IfL + QfR QfR - m f(qfRQfL + qfLQfR)} (127) 

/=i 

qjL and qfR can in general not be rotated separately any longer because of the quark masses. The actual 
flavour symmetry therefore depends on the quark mass matrix 

M q = diag (m u , m d , m s , m c , m b , m t ) . (128) 

In order to find all even only approximate symmetries of QCD, we distinguish several cases de- 
pending on the specific values of the quark masses. 

i. In the real world, all quark masses are non-zero and they are all different from each other. In this 
case, the remaining flavour symmetry amounts to the phase transformations qf lr — ► e~ l£ f qj l^r 
(f = 1, . . . , 6) where the phase Sf for a given flavour must be the same for qji and qfR. The 
symmetry group reduces to the product U{\) x U(l) x ••• x U(l) = U(l) 6 leading to the well- 
known conserved flavour quantum numbers N u , Nd, N s , N c , iVj, and N t . All these symmetries are 
broken by the weak interactions, except their sum (baryon number) 

B = (N u + N d + N s + N c + N b + N t ) /3 . (129) 
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ii. In some approximation, the quark masses are still non-zero but uf of them are equal (np < Np = 
6). In this case, the maximal symmetry group U(6) x U(6) reduces to 

U{n F ) x U(lf- nF ~ SU(n F ) x U(l) x U(lf- nF . (130) 

The only realistic cases are np = 2 or 3 and they lead to well-known approximate symmetries: 

np = 2: m u = ma — ► isospin SU(2) 

np = 3: m u = m^ = m s — ► flavour SU(3) 

iii. A much more radical approximation consists in setting some of the quark masses to zero: mj = 
(/ = 1, . . . , np). In this case, np of the qji and q^ can again be rotated separately implying the 
chiral symmetry 

SU(n F ) L x SU(n F ) R x U(l) v x U(1) A [xU(lf~ nF ] . (131) 

To set np = 2 of the quark masses to zero is a reasonable approximation in view of m U: d *C 
Aqcd» whereas np = 3 (setting also m s = 0) is certainly more daring. U(l)v is again responsible 
for baryon number conservation. The factor U(1)a is actually not a symmetry of full QCD at the 
quantum level (abelian anomaly). 

We are familiar with isospin and flavour 577(3) that we see at least approximately realized in the hadron 
spectrum. But what are the consequences of the approximate chiral symmetry of QCD? If chiral symme- 
try would manifest itself in the hadron spectrum, each hadron would have to have a partner of opposite 
parity of approximately the same mass. That is obviously not the case: chiral symmetry appears to be 
more hidden than isospin, for example. In order to understand the manifestations of chiral symmetry, we 
have to recall the main features of 



Spontaneous symmetry breaking 



There are many examples of SSB in physics and the best-known example in particle physics is the 
spontaneously broken electroweak symmetry (see the Lectures of W. Buchmuller at this School). 

The mechanism was first realized in condensed matter physics and a good example for our purpose 
is the ferromagnet. The underlying theory of the ferromagnet (eventually QED) is certainly rotationally 
invariant. The Hamiltonian (e.g., of the Heisenberg model) does not single out any direction in space. 
Nevertheless, in the ground state of the ferromagnet the spins align in a certain direction. The direction 
is arbitrary and there is no trace of it in the Hamiltonian. In this sense rotational symmetry is "sponta- 
neously" broken. It is certainly not manifest in the ground state but it has other important consequences 
such as the existence of excitations (called magnons or spin waves in this case) with a very special disper- 
sion law. The dispersion law is the dependence of the frequency uj on the wave length A or on the wave 
number k = 2tt/\. SSB in condensed matter physics implies that for some excitations the frequency 
tends to zero for infinite wave length: 

limw(fe) = lim u(k) = . (132) 

This property of magnons in particular is easy to visualize as shown in Fig. 23. In the ground state 
of the ferromagnet spins are aligned. A typical spin wave is displayed in the second line: the wave length 
is the distance between spins pointing in the same direction. In the limit A —> oo, the spins become again 
aligned, albeit in a different direction in general. Since by the assumed rotational symmetry of the theory 
each direction is as good as any other, the configurations in the first and in the last line must have the 
same energy as expressed by Eq. (132). Put in another way, magnons do not have an energy gap in their 
spectrum. 
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Fig. 23: Spin directions in the (one-dimensional) ferromagnet: ground state, spin wave with finite wave length and 
a spin wave with infinite wave length. 

What is the analogy in a relativistic QFT like QCD? The frequency is replaced by the energy of the 
particle, with the three-momentum instead of the wave number. The dispersion law is just the relativistic 
energy-momentum relation E = \Jp 2 + m 2 . The energy tends to zero for p — ► if and only if the 
particle is massless: 

lim E = lim Vp 2 + m 2 = <^=^ m = . (133) 

How can we distinguish if a symmetry is realized in the usual way (like isospin) or if it is spontaneously 
broken? The crucial question is what the conserved charge Q = f d 3 xJ°(x) associated with a symmetry 
current J fl (x) (d^J^ = 0) does when applied to the ground state (vacuum). In a relativistic quantum 
field theory, there are only two options: 

Goldstone alternative 

Q|0) = ||Q|0)|| = oo 

Wigner-Weyl Nambu-Goldstone 

linear representation nonlinear realization 

degenerate multiplets massless Goldstone bosons 

exact symmetry spontaneously broken symmetry 

The left column describes the more familiar case (Wigner-Weyl) where states are grouped in multiplets 
(irreducible representations of the symmetry group). The vacuum is annihilated by the charge and the 
states in a given multiplet all have the same mass (degeneracy). In the other possible realization (Nambu- 
Goldstone), applying the charge to the vacuum is, strictly speaking, not defined. There are no degenerate 
multiplets (therefore we don't see the symmetry in the spectrum) but there must be massless particles 
in the theory (Goldstone bosons). Although the charge cannot be applied to the vacuum directly, the 
following matrix element called an order parameter of SSB may be well defined: 

<0|[Q,A]|0) (134) 

where A is some operator. If we can find an order parameter that is different from zero, the symmetry 
associated with the charge Q is necessarily spontaneously broken. This is easy to understand: because 
of the commutator in (134) the order parameter vanishes if Q annihilates the vacuum. In the familiar 
example of electroweak symmetry breaking, scalar field operators ipi take the place of A with [Q, tpi\ = 
Cijifj. If Cij(0\ifj\0) / (Higgs vacuum expectation value), the electroweak symmetry is spontaneously 
broken. This is also a good example that the argument goes in one direction only. Even if all scalars 
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have zero vacuum expectation values or if there aren't any scalar fields at all the symmetry may still 
be spontaneously broken. The mechanism could involve some other operator A in the order parameter 
(134). 

For each spontaneously broken symmetry Goldstone's theorem implies the existence of a massless 
state |G) with 

(0|J°(0)|G)(G|A|0) /0. (135) 

A necessary and sufficient condition for SSB is that the 



Goldstone matrix element (0|J°(0)|G) / 



implying also that the Goldstone state |G) has the same quantum numbers as J°(0)|0). The following 
remarks are useful: 

• The state | G) need not correspond to a physical particle. This can only happen in the case of a 
spontaneously broken gauge symmetry as in electroweak theory. 

• J°(0) is usually a rotationally invariant bosonic operator and thus |G) carries spin 0. 

• Spontaneous breaking of discrete symmetries does not give rise to Goldstone bosons. 

The main features of SSB can be discussed in the original Goldstone model. It has a single complex 
scalar field tfi{x) with the Lagrangian 

v 2 
^Goldstone = d^d^cfy - A(00 t - — ) 2 (A, v real and positive) , (136) 



v 2 



symmetric with respect to global U(l) transformations (f>(x) — ► e lOL 4>{x). The minimum of the Mexi- 
can hat potential 
R(x),G(x) with 



can hat potential occurs at ipcfr = — . Decomposing the complex field (j)(x) into two hermitian fields 



<f>{x) = (R{x)+iG{x))/V2 (137) 

(0\R(x)\0)=v, (0|G(x)|0) =0, 

the Lagrangian expressed in terms of the fields R(x), G{x) displays the following spectrum at tree level: 

Goldstone boson field G{x) M G = 

massive field H (x) = R(x) — v Mh = \/2\v 

Denoting the four-momenta of Goldstone particles generically as pq, one finds an unexpected behaviour 
for scattering amplitudes: they vanish for p G — ► 0, e.g., 

A(GG -► GG) = 0(p%), A(GH ->■ GH) = 0(p 2 G ) (138) 

for arbitrary values of non-Goldstone momenta pn- More generally, Goldstone bosons decouple when 
their energies tend to zero. 

This behaviour looks mysterious at first, but it can be understood by choosing a different set of 
fields. Instead of the fields G(x) and H(x), we choose another representation (polar decomposition) that 
may be familiar from electroweak theory: 

<f>( x ) = -L[/i(x) + v]e i 9( x )/ v . (139) 

v2 
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In terms of the hermitian fields g(x), h(x) the Goldstone Lagrangian takes the form 



\^9? + ^ 2 



Aioidstonc = ^g) 2 + -^(h 2 + 2vh)(d IM g) 2 (140) 



+ -(d lt h) 2 -Xv 2 h 2 --{h A + 4vh 3 ). 

A general theorem of QFT ensures that the fields G, H on one side and g, h on the other side produce 
the same S'-matrix elements although the Green functions are in general different. 

The main consequences are the following. 

- The particle spectrum is unchanged: 

Goldstone field g{x) M g = 

massive field h(x) Mh = \J2\v 

- The Goldstone field g has only derivative couplings implying for the scattering amplitudes consid- 
ered previously: 

lim A(p G ) = . (141) 

The important lesson is very general and not restricted to the Goldstone model. S'-matrix elements 
with only Goldstone states vanish for pa — ► 0. When other non-Goldstone particles participate in the 
initial and final states, the statement remains true for some matrix elements like for elastic scattering 
G H — ► G H. In general, the interactions of Goldstone bosons among themselves and with other matter 
become arbitrarily weak for small momenta. 

3.5 Chiral perturbation theory 

We start from a theorist's paradise (copyright H. Leutwyler), QCD in the chiral limit where np = 2 [or 
3] quarks u, d [, s] are massless: 

^QCD = QL0QL + qR0QR + -Cheavy quarks + £ga U ge ( 142 ) 

with 

q T = (ud [s]) . 

As explained in the previous section, this Lagrangian has a global symmetry 

SU{n F ) L x SU(n F ) R x U(l) v x 17(1)^ [xU(lf- nF ] . (143) 

The nonabelian factor G = SU(riF)L x SU{tif)r is called the chiral group. 

Although not yet proven from QCD alone, there is strong evidence, both from phenomenology and 
from theory, that chiral symmetry is spontaneously broken. The spontaneous breaking does not affect all 
of G but, roughly speaking, half of it: G — ► H = SU{np)v- The so-called vectorial subgroup H of 
G is nothing but isospin (for rip = 2) or flavour 577(3) (for np = 3) and it is realized in the familiar 
way a la Wigner-Weyl. Some arguments in favour of this spontaneous breakdown are: 

As already emphasized before, there are no parity doublets in the hadron spectrum. 

• There is no other convincing argument why the pseudoscalar mesons are the lightest hadrons. 
Spontaneous chiral symmetry breaking implies that they would be massless in the chiral limit 
(pseudo-Goldstone bosons). 

• The vector and axial-vector spectral functions are very different (p vs. a{). 

• The so-called anomaly matching conditions together with confinement require that G is sponta- 
neously broken for np > 3. 



• Under very reasonable assumptions, SU(riF)v is not spontaneously broken. It is of course expli- 
citly broken by the differences between quark masses. 

Even if it has not been possible so far to prove directly from QCD that chiral symmetry is sponta- 
neously broken, we can ask for possible order parameters. It turns out (more details can be found 
in Ref. [33], for instance) that the simplest such order parameter involves the pseudoscalar operators 
Ab = q^Xbq {a = 1, . . . ,8) giving rise to the quark condensate (0\qq\0). There is evidence both from 
lattice QCD and from phenomenology that the quark condensate is non-zero implying spontaneous chiral 
symmetry breaking. As will be discussed in Sec. 3.7, the quark condensate is in fact the dominant order 
parameter of spontaneous chiral symmetry breaking, in a sense to be specified later. 

From Goldstone's theorem we know (still in the chiral limit) that there are np — 1 massless Gold- 
stone bosons: 

np np — 1 Goldstone bosons 



2 3 7T 

3 8 n,K,n 



Although the real world is not a theorist's paradise, we still expect low-energy amplitudes to be domi- 
nated by the exchange of pseudoscalar mesons, which are the lightest hadrons also in the real world. In 
order to calculate such amplitudes, we construct an effective field theory with only Goldstone fields. As 
already explained in Sec. 3.1, the Lagrangian of Goldstone fields is nonrenormalizable and it is in fact 
even nonpolynomial. The underlying physical reason is that we can add any number of sufficiently soft 
pions (still massless!) to a hadron state without appreciably changing its energy. Therefore, we have 
degenerate states with different numbers of particles that are related by chiral symmetry transformations. 
For the Lagrangian this argument implies that the symmetry transformations are nonlinear in the pion 
fields. Starting from a Lagrangian with fixed powers in the Goldstone fields, successive nonlinear trans- 
formations generate any number of fields in the Lagrangian. Since the Lagrangian is to be symmetric 
under such transformations it must necessarily be nonpolynomial. 

The basic building block of chiral Lagrangians is therefore a nonpolynomial matrix function of 
the Goldstone fields, e.g., the exponential function (for np = 3) 



U(<f>) = exp(iV2$/F), $ 
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where F is the pion decay constant in the chiral limit characterizing the size of the Goldstone matrix 
element (0|J°(0)|G). 

Chiral Lagrangians are organized according to the number of derivatives of the fields. The unique 
lowest-order Lagrangian of 0(p 2 ) with two derivatives is the so-called nonlinear a model: 

4 0) = ^n F (dpUdPlfl) =: ^-(d^Ud^) = d^+d^ir- + \d^d^ + 0(vr 4 ) , (145) 

using a bracket notation for np -dimensional traces. 

So much for the paradise. Back to reality, we must admit that there is no chiral symmetry in 
nature! In the Standard Model, it is explicitly broken in two different ways. 

• Chiral symmetry is explicitly broken by nonvanishing quark masses. This should be a small mod- 
ification for two, a more pronounced one for three flavours: 
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m u , md *C M p rip = 2 

m s < M p tif = 3 

• Also the electroweak interactions break chiral symmetry. If electroweak effects are to be included, 
they can be taken into account perturbatively in a, Gf- 

The main assumption of chiral perturbation theory (CHPT) is that an expansion around the chiral 
limit (the theorist's paradise) makes sense. Therefore, even in the absence of electroweak interactions, 
chiral Lagrangians are organized in a two-fold expansion. 

i. Spontaneous chiral symmetry breaking gives rise to an expansion in derivatives of the fields leading 

to an expansion of amplitudes in the momenta of pseudo-Goldstone bosons, 
ii. Explicit symmetry breaking suggests an expansion also in the quark masses m q . 

The two expansions can be related via the meson masses. As will be discussed in the next subsection, 
the squares of the meson masses start out linear in the quark masses: 

M 2 M ~m q + 0{m 2 q ) . (146) 

The standard chiral counting therefore amounts to treating quark masses like the second power of mo- 
menta: 

m q = 0(M 2 M ) = 0{p 2 ) . (147) 

The effective Lagrangian (for pseudoscalar mesons) is therefore of the form [34] 



£eff — C2 + C4 + Cq + . . . 

F 2 
4 



£2 = —(dpUd^+xU^ + X^U) (148) 



where \ represents the quark masses: x = ^BM. q = IB diag(m u ,md[, m s ]). The lowest-order La- 
grangian contains only two parameters F, B that are related to physical quantities as 

F n = F [1 + 0(m q )} , (0\uu\0) = -F 2 B [1 + 0(m q )] . (149) 

The lowest-order amplitudes of CHPT are of 0{p 2 ) and they correspond to the current algebra 

amplitudes of 40 years ago. The tree-level amplitudes can be read off directly from the Lagrangian (148) 

depending only on F n and M^ [M 2 = B(m u + m<f), •••]• For instance, the elastic tttt scattering 

amplitude of 0(p 2 ) is given by 

s — M 2 
A 2 (s,t,u) = —^. (150) 

* IT 

Contrary to symmetries like isospin that relate different amplitudes, the spontaneously broken chiral 
symmetry makes an absolute prediction for this scattering amplitude. It was left as an exercise to the 
audience to explain why that is possible. 

The lowest-order results we have been discussing so far were already known in the late sixties and 
early seventies of the last century (current algebra, phenomenological Lagrangians). After an influential 
paper of Weinberg [35], but especially with the work of Gasser and Leutwyler [34] the systematic treat- 
ment of QCD at low energies became a respectable theory. The first step was to construct the Lagrangian 
of next-to-leading order £4 that contains 10 (7) additional coupling constants (usually called LECs for 
low-energy constants) for SU(3) (SU(2)). With a hermitian Lagrangian tree amplitudes are necessarily 
real but we know that unitarity and analyticity require complex amplitudes. It is not difficult to convince 
oneself that imaginary parts occur first at 0(p 4 ). The consequence is that a systematic low-energy ex- 
pansion requires a loop expansion beyond lowest order [35]. But loop amplitudes have a tendency to be 
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divergent. For a bona fide QFT we therefore need both regularization and renormalization. As strange 
as it may sound, also nonrenormalizable theories can and in fact must be renormalized to qualify as 
respectable QFTs. 

A nonrenormalizable QFT like CHPT has many common features with the more standard renor- 
malizable theories. 

•* Divergences are absorbed by the coupling constants in the higher-order Lagrangians C\,C§,.... 
Unlike in renormalizable theories, new LECs occur at every order of the chiral expansion. 

•* The renormalized LECs are scale dependent just like the strong coupling constant g s (n). They 
can be interpreted as describing the effect of all heavy hadronic states that are not represented by 
explicit fields in the Lagrangian. 

<•" Renormalization ensures that there is no dependence on some artificial cutoff. 

For phenomenological applications, we have to know the values of the various LECs. In principle, QCD 
fixes those constants but a matching between QCD and CHPT is not possible in perturbation theory. 
This was already discussed in Sec. 3.1 for general EFTs with SSB but it is also easy to understand in the 
present case: CHPT can only be applied for energies E <C M p whereas perturbative QCD only makes 
sense for E > 1.5 GeV. The most successful approaches bridging this gap to get information on the 
LECs are resonance saturation (based on the properties of QCD for large N c , i.e. for a fictitious world 
with many colours) and lattice QCD. 

The chiral expansion is an expansion in p 2 /(AttF) 2 where p is a characteristic momentum for the 
process in question. Therefore, the chiral expansion should and does work better for SU{2) than for 
SU{3): 

»'- 2: (^i- MM ^. B '- S! (iSF- au, jJijr- (151) 

Most amplitudes and form factors for realistic processes have been calculated at least to next-to- 
leading order. There is an easy-to-use Mathematica program to generate both strong and nonleptonic 
weak transitions up to 0(p A ) that is described in Ref. [36] and is available for general use. The state of 
the art is next-to-next-to-leading order or 0(p 6 ). A short introduction can be found in Ref. [33]. 

3.6 Light quark masses 

In CHPT, the light quark masses always appear in the combination B m q ~ m q (0|u u|0). Since there 
are no quarks or gluons in CHPT, only QCD scale invariant quantities can appear. As we discussed in 
Sec. 3.2, quark masses are scale dependent whereas the product m q (0|uu|0) is not. The consequence is 
that CHPT can only provide methods for extracting the ratios of quark masses. 

The lowest-order expressions for the meson masses in terms of quark masses can be read off 
directly from the Lagrangian (148): 

M 2 + = 2mB , M 2 = 2mB + O [{m u - m d ) 2 /{m s - m)] 

2 
M 2 K+ = (m u + m s )B , M 2 S = -{m + 2m s )B + O [(m„ - m d ) 2 /{m s - m)] 

M^o = (m d + m s )B , m := -{m u + m d ) ■ (152) 

Several well-known relations follow from these expressions: 



40 



Gell-Mann-Oakes-Renner 
GMOR, Weinberg 
Gell-Mann-Okubo 



FlMl 



-2m(0\uu\0) 



B 



Ml 

2rh 



M K + 



M 2 K0 



m s + m u m s + m d 

(isospin limit) 

The relations (152) are also the basis for the so-called current algebra mass ratios 
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These ratios are subject to higher-order corrections, most importantly of 0(p 4 ) = 0{m? q ) and 0(e 2 m s ). 
Because of an accidental symmetry at 0(p A ), the ratios to s /to^, m u /md cannot be extracted separately 
from S-matrix elements but only in the combination known as Leutwyler's ellipse [37] shown in Fig. 24. 




Fig. 24: Constraints on light quark mass ratios [37]. 



In addition to the full boundaries following directly from CHPT (the difference is due to uncertainties in 
the electromagnetic corrections), information is also available from 77 — 77' mixing (dotted lines), baryon 
mass splittings and p — lo mixing (dashed lines). The overall conclusion is that the corrections of 0(p 4 ) 
are small for the ratios. The next-to-next-to-leading corrections of 0(p 6 ) are also known [38] but there 
are at the moment too many unknown LECs for quantitative conclusions. 



0{p 2 ) 

o( P 4 ) 



m u /m d 



m s /m d 



m s /rh 



0.55 
0.55 ± 0.04 



20.1 
18.9 ±0.8 



25.9 
24.4 ± 1.5 



Table 1: Quark mass ratios to 0(p 4 ) [37]. 

Absolute values of the quark masses are less well known than the ratios. The main methods are QCD sum 
rules and lattice simulations, most recently with full (unquenched) QCD. From the Review of Particle 
Properties [14], the combined result of lattice determinations of the strange quark mass is m s (2 GeV) = 
(100 ± 25) MeV. More results are reproduced in Fig. 18. 
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3.7 Pion pion scattering 

Pion pion scattering has a privileged status in CHPT. It is the fundamental scattering process of CHPT 
and it involves only pions. The low-energy expansion can therefore be set up in the framework of 
chiral 51/(2) and it can be expected to converge well. The scattering amplitude is very sensitive to the 
mechanism of spontaneous chiral symmetry breaking giving access to the quark condensate in particular. 

The following review of recent developments will be restricted to the isospin limit (m u = md) 
in the absence of electromagnetic corrections. In this case, the information for all possible scattering 
channels is contained in a single amplitude A(s, t, u) (with s + t + u = 4M^). 

The lowest-order amplitude of 0(p 2 ) was already shown in Eq. (150): 

a r , n S ~ M - 

A 2 {s,t,u) -- 



F 2 ' 

x IT 

m 1M 2 
At the same order, the quark mass ratio r = — - = if — 1 ~ 26, as also shown in Table 1. As the 

M rh Ml 

mass ratio r, also the S-waves are very sensitive to the quark condensate. In a modified version of CHPT 
(Generalized CHPT [39]), one can tune the quark condensate. As an example, I show the leading-order 
results for the 1 = S-wave scattering length a[] and for the quark mass ratio r for both the standard and 
for a very small value of the quark condensate: 

qg r B{y = 1 GeV) 

016 26 1.4 GeV (standard value) 
026 10 Fk 

At next-to-leading order, the scattering amplitude was calculated in 1983 [40]. It turns out that especially 
the S-wave scattering lengths are quite sensitive to chiral corrections (chiral logs). For instance, a[] 
increases from 0.16 at 0(p 2 ) to 0.20 at 0(p 4 ), an increase of 25 % and thus quite a bit larger than the 
natural estimate in Eq. (151). Since the favoured experimental value of a(j at that time was 0.26 (with a 
25 % error), it seemed mandatory to perform one more step in the chiral expansion. From the amplitude 
to 0(p 6 ) [41] it was clear that a value aj] = 0.26 was not compatible with QCD. To narrow down 
the uncertainties related to the LECs appearing in the -kit amplitude, the chiral amplitude was finally 
combined with dispersion theory (Roy equations) [42]. 

CHPT together with dispersion theory predicts not only the S-wave scattering lengths with amaz- 
ing precision [42], 

ag = 0.220 ± 0.005 , ag = -0.0444 ± 0.0010 , (154) 

but also the S- and P-wave phase shifts. As a by-product, the pionium lifetime is predicted to be r = 
(2.9 ±0.1) • 10~ 15 . A short history of the S-wave scattering lengths is shown in Fig. 25. 

There is a small caveat here. All the results have been derived in the standard framework of CHPT 
that assumes a large quark condensate. With recent experimental information on the pion-pion phase 
shift difference 5® — 5{ from K e 4 decays, even this last loophole could be closed. Using the correlation 
between Oq and a(] implied by the Roy equations, the measured phase shift difference [43] can be used 
to determine a(] as shown in Fig. 26. The fitted value [44] a[j = 0.221 ± 0.026 is in perfect agreement 
with the prediction (154) from the combined analysis of CHPT and Roy equations. 

The precise determination of the tttt scattering amplitude has several important implications. One 
application concerns the chiral expansion of the pion mass in terms of the light quark masses m u , rrid'. 

Ml = M 2 --^ IW1 M^ + 0{M % ) (155) 

M 2 = {m u + m d )\(0\uu\0)\/F 2 . 
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A by-product of the analysis of mr scattering is a precise determination of LECs like £3, which implies in 
turn that more than 94 % of M w are in fact due to the leading term (the Gell-Mann-Oakes-Renner term) 
confirming the standard mechanism of spontaneous chiral symmetry breaking [44]. In other words, the 
quark condensate is indeed the dominant order parameter of chiral symmetry breaking. 
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Fig. 25: Predictions for S-wave scattering lengths 
from current algebra till today, taken from Ref. [42]. 
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Fig. 26: K e 4 data translated into a determination of 
the 1 = S-wave scattering length a§ [44]. 



The precise knowledge of the -kit scattering amplitude from CHPT and dispersion theory has more 
recently (after the School) produced another important insight. The much discussed scalar isoscalar state 
now called /n(600) by the Particle Data Group [14], but more commonly known as the a meson, was 
analysed by Caprini, Colangelo and Leutwyler [45] by extending the Roy equations to complex values 
of the Mandelstam variable s. The result is an astonishingly precise determination of mass and width of 
the lightest hadronic resonance: 
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The a resonance has the quantum numbers of the vacuum and it corresponds to an unambiguous pole on 
the second sheet of the scalar isoscalar partial wave. Its real part is close to threshold but the pole is quite 
far from the real axis. Although not as straightforward as its position in the complex energy plane, the 
most appealing interpretation of the a is a quasi-bound state of two pions, quite different from a member 
of a standard q q meson nonet [42, 45]. 



3.8 K13 decays and V us 

The Cabibbo-Kobayashi-Maskawa (CKM) matrix Vij determines the structure of the hadronic charged 
weak current. The matrix elements are fundamental parameters of the Standard Model. Together with the 
quark and lepton masses, the CKM matrix and the corresponding lepton mixing matrix contain informa- 
tion about the mechanism of mass generation. Anticipating forthcoming LHC data to unveil the secrets 
of mass generation (Higgs sector), both masses and mixing angles should be measured as precisely as 
possible. 

With three generations of quarks, the CKM matrix must be a unitary matrix. For some time, a 
possible problem with CKM unitarity has been discussed. With the PDG values of 2004 [14], 

\V ud \ = 0.9738(5) , \V US \ = 0.2200(26) , (157) 

unitarity appeared to be violated at the 2.2 a level by the elements of the first row V u j (j = d,s,b): 

1 2 
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V, 
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-0.0033(15) . 



(158) 
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At this level of accuracy, the element V u b is still negligible. On the other hand, the required precision for 
Vud and V us calls for reliable isospin violating and electromagnetic corrections. 

The most accurate determination of V us , both experimentally and theoretically, comes from Ki% 
decays that can be treated in the framework of CHPT. The decay amplitude is governed by two form 
factors f+(t) and f-(t) with t = (px — P-k) 2 ' 

{<k-(Pk)\siA K °(pk)) = f+ >7T ~(t) (PK+Pir)^ + /f ° w "(t) (PK-Pn)n ■ (159) 

Both form factors are known to next-to-next-to-leading order in CHPT. For the extraction of V us , the 
crucial quantity is /+(0). The chiral expansion is of the form 

/*°*~ (0) = 1 + / p4 + f e2p 2 + f p e + 0[(m u - m d )p\ e 2 p 4 } . (160) 

The status of the various contributions is as follows: 

fpi —0.0227 (no uncertainty) Gasser, Leutwyler [46] 

/ e 2 p 2 radiative corrections Cirigliano, Neufeld, Pichl [47] 

f p 6 loop contributions Bijnens, Talavera [48]; Post, Schilcher [49] 

tree contributions LECs L\, Cu + C34 



As a first test, we look at the ratio 



1/2 



r+0 \nK m )Mi + i K+ ) \ff--m' { } 

This ratio is independent of f p a in Eq. (160) and it can therefore be predicted quite accurately [47,50]: 

r% = 1.023 ± 0.003 , (162) 

to be compared with the experimental value [51]: 

r+P = 1.036 ± 0.008 . (163) 

What could be the origin of a possible discrepancy that is also suggested by the compilation of recent 
data in Fig. 27? 

• In the past, radiative corrections have not always been state of the art. Nowadays, also experimen- 
talists should only use the modern CHPT treatment [47]. 

• Measurements of the K + and Kl lifetimes should still be improved. 

• On the theory side, an unlikely but in principle still possible explanation could be that the error 

due to effects of 0[(m u — m^p 4 , e 2 p 4 ] is underestimated. 

The contribution /„6 is the sum of a loop and of a tree-level part: 
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= 0.0093 ± 0.0005 


Bijnens, Talavera [4( 
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Fig. 27: Present experimental and theoretical status of /+(0) \V US | taken from Ref. [53]. 



The last equation is based on a large- N c estimate of Ref. [52]. As can be seen in Fig. 28, two terms 
interfere destructively weakening the overall dependence on the scalar resonance mass Mg. The same 
interference leads to a rather modest scale dependence of the result for M v < \x < 1 GeV. 




M P = 1.3 GeV 



1.3 



1.5 



M s (GeV) 
Fig. 28: Tree-level contribution of 0(p 6 ) to f? " (0) and its two components [52]. 



The final result, including the uncertainty due to a possible second pseudoscalar multiplet P' , 



is [52] 



fT(M p ) 



f 



p> 
p 6 

(0) 



fv* 



-0.002 ± 0.008 
0.007 ±0.012 
0.984 ±0.012 . 



l/N< 



_i_ n nno +0.000 

±U.U02 M s -0.002 P' 



(166) 



With our large- N c estimate of the tree-level contribution of 0(p 6 ), f p a is dominated by the loop part 
(164). It exhibits less SU(3) breaking than the well-known result of Leutwyler and Roos [54] and 
a recent lattice estimate [55]. Taking the most recent value of V u d and assuming unitarity of the CKM 
matrix, the predictions can be compared directly with the experimental results as shown in Fig. 27. A new 
result for the neutron lifetime [56] would imply a shift of all theoretical values in Fig. 27 to the left but 
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the corresponding accuracy of V U( i is not yet competitive with super-allowed nuclear Fermi transitions. 
Another way to read Fig. 27 is that the estimate of Ref. [52] leads to a value of V us , 

\V US \= 0.2208 ± 0.0027 /+(0) ± 0.0008 cxp , (167) 

that is a bit smaller than the unitarity prediction. 

Finally, the slope of the scalar form factor that also depends on the LECs C12, C34 can also be 
predicted [52] and it is in good agreement with the recent most precise experimental value [57] : 

A = (13 ± 3) • 10~ 3 Cirigliano et al. [52] 

(168) 
A = (13.72 ± 1.31) • 10~ 3 KTeV [57] 

4 Summary and epilogue 

There is an amazing richness contained in the simple Lagrangian that we "derived" from the existence 
of colour and from the gauge principle: 

N F 



£qcd = -\tr(G^G» u ) + J^q f [0 - m f l c ) q f 



There is in fact no better summary for these lectures. 

It may seem a long way from the naive quark model to QCD but it all happened in less than ten 
years. On the asymptotically free side, perturbative QCD is a complete success and it will be especially 
needed for understanding the background for new physics at the LHC and beyond. 

There is much more left to be understood at the other end of the scale. If confinement is really 
forever, we would witness the first case in the history of physics when new constituents of matter have 
been identified beyond reasonable doubt and yet they can never be isolated. The question sounds prepos- 
terous but it is well supported: have we already arrived at the innermost structure of hadrons? Or put in a 
different way, is there no further structure to be expected before strings and quantum gravity eventually 
take over? In a short time, once the LHC will produce first results, we may learn how to rephrase the 
question. 
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